NNLO Logarithmic Expansions and Precise Determinations of 

the Neutral Currents near the Z Resonance at the LHC: 

The Drell-Yan case 



O 

O . 3 

(N 



^'^Alessandro Cafarella, ^'^Claudio Coriano, ^Marco Guzzi 

^Institute of Nuclear Physics, NCSR "Demokritos" , 15310 Athens, Greece 

Dipartimento di Fisica, Universitd del Salento 
and INFN Sezione di Lecce, Via per Arnesano, 73100 Lecce, Italy 



Department of Physics and Institute of Plasma Physics, University of 
Crete, GR-710 03 Heraklion, Crete, Greece 



'sj- ■ Abstract 
> " 

'. We present a comparative study of the invariant mass and rapidity distributions 

, in Drell-Yan lepton pair production, with particular emphasis on the role played 

O . by the QCD evolution. We focus our study around the Z resonance (50 < Q < 

. 200 GeV) and perform a general analysis of the factorization/renormalization scale 

dependence of the cross sections, with the two scales included both in the evolution 
and in the hard scatterings. We also present the variations of the cross sections due 
Oh! to the errors on the parton distributions (pdf 's) and an analysis of the corresponding 

. X-factors. Predictions from several sets of pdf's, evolved by MRST and Alekhin 

are compared with those generated using Candia, a NNLO evolution program 
^ ■ that implements the theory of the logarithmic expansions, developed in a previous 

■ work. These expansions allow to select truncated solutions of varying accuracy using 

the method of the x-space iterates. The evolved parton distributions are in good 
agreement with other approaches. The study can be generalized for high precision 
searches of extra neutral gauge interactions at the LHC. 
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1 Introduction 



Accurate determinations of the QCD background at the LHC, especially for some selected 
hadronic cross sections, are going to be very important in order to increase our potential for 
new discoveries. For this reason it is necessary to know the size of the radiative corrections 
to some selected processes at higher orders. At the same time, the quantification of the 
impact of the errors in the determination of these observables is going to be critical 
in order to enhance our confidence on the reliability of the perturbative expansion. It 
is particularly so in the search for extra Z', which are ubiquitous in extensions of the 
Standard Model [2] for instance in models derived from the string construction 
[Ij or with extra neutral interactions modified by an anomaly inflow [5], where mass- 
dependent and anomaly-related corrections require very high accuracy to be properly 
identified and separated from the large QCD background. 

With these motivations in mind we have proceeded with an independent analysis of 
the next-to-next-to-leading order (NNLO) QCD corrections, starting from an accurate 
investigation of the impact of the evolution on the physical observables at the LHC. In 
this context, the role played by the Drell Yan cross section is particularly important. 
In fact, the possibility of discovering extra neutral currents at the new collider may be 
related to the determination of this cross section with very high accuracy far beyond 
the Z peak, at values of the invariant mass of the lepton pair up to 5 TeV [T], which 
is commonly thought to be the upper limit in searches of this type. For this reason, a 
precise determination of the pdf 's at any value of the Bjorken variable a; [6j is needed (see 
[7j and refs. therein). Issues of resummation of the perturbative expansions become also 
critical for the correct determinations of several distributions at the edge of phase space 
(see for example [H]). 

At this time, while our knowledge of the role played by the coefficients of the QCD 
hard scatterings in some key partonic processes is quite satisfactory, that of the behav- 
ior of the pdf 's is not of a comparable level, and the model-dependence of the various 
parameterizations is still large. Recent parameterizations of the pdf 's come with the quan- 
tification of their errors, presented to next-to-leading order (NLO) and NNLO, whenever 
possible, obtained by the fits used by various groups (we limit our analysis to [9J and 
[TU] ) to match several sets of experimental data in pp collisions, such as Drell- Yan and 
Deep Inelastic Scattering (DIS). These errors, which estimate the goodness of a fit, are 
naturally thought of being of experimental origin. But there are also other sources of 
indetermination, mostly of theoretical origin, which need to be taken into consideration. 
One of them, apparently of more trivial nature, is related to the way the evolution is 
implemented through NNLO. 
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At a first look this last point might be misinterpreted and the corresponding "error" 
coming to be attributed to the "model dependence" of a given parameterization set, while 
it amounts to a theoretical indetermination, intrinsic to perturbation theory, since it is 
going to be there for any chosen model of pdf's. The reason is simple and also quite 
immediate: there is not a unique approach to solve the DGLAP equation, and, again, 
not for a numerical/algorithmic reason, related to the limited numerical precision of a 
given algorithm. In fact, a given solution, of a typical accuracy, organizes the logarithmic 
corrections in a specific and unique way. These are summed or, eventually, resummed 
if exact or, instead, accurate (truncated) solutions are selected. Therefore, the issue of 
determining the best possible way to solve the evolution does not seem to have a unique 
answer, being directly related to the possibility of choosing among different theoretical 
approaches, all equally acceptable. The goal of this work is threefold: to test the accuracy 
of the logarithmic expansions proposed by us in a previous work by comparing with other 
methods of solution; to quantify the experimental errors on the DY cross section coming 
from the pdf's and to analyze the differences between accurate and "exact" solutions of 
the evolution equations. The method used by us to solve the DGLAP equations is based 
on a simple perturbative organization of the logarithmic corrections that we call "the 
NNLO logarithmic expansion in x-space" . These expansions work for kernels given to all 
orders. 

Thanks to the availability of the two-loop evolution kernels [H] and of the correspond- 
ing hard scattering coefficients, which in the Drell-Yan case have been known for some 
time [12], complete NNLO analysis are now possible and allow to perform sophisticated 
tests that give to us the opportunity to quantify the effects that we have just mentioned. 
In the case of Drell-Yan both the total cross section and, with a modification, the rapid- 
ity distributions of the lepton pair on the final state are at reach. The latter have been 
presented recently [13\, together with a dependence of the predictions on the factoriza- 
tion scale, which is important in order to monitor the overall stability of the perturbative 
series in a^, the strong coupling constant. This requires the determination of the cross 
section for a varying factorization scale and of the relevant i^"-factors at various energies. 
However, in order to determine in a robust way the size of the NNLO corrections and 
the role that they play in some important predictions, we believe that it is mandatory to 
perform an independent analysis of the evolution, defining benchmarks for the evolution 
of the pdf's from different perspectives respect to the common ones. These are based 
either on numerical discretizations (so called brute force methods), which are affected by 
contributions of all-orders in the strong coupling, or on the methods of the Mellin mo- 
ments, using special expansions. Accuracy means that we define the perturbative solution 
so to include only parts of the corrections, in a certain expansion in the strong coupling 
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|14] . given our limited knowledge of the perturbative expansion of the kernels. 

It is well known that the issue of accuracy in the choice of the solution of the equation 
has never been fully addressed in the previous literature. While this issue is less important 
at NLO, given the size of the i^-f actors which are about 20% in the region that we explore 
and at the energy that we select (v^ = 14 TeV) for a p-p collider, things become more 
subtle at NNLO, where the relative ii"-factors relating the NLO to the NNLO cross section 
are determined to be much smaller. We will show that it is the QCD evolution to drive 
the NNLO cross section to an overall reduction in the region that we have analyzed. As a 
result of our analysis we are able to quantify the theoretical error implicit in the various 
choices of the evolution scheme, which are smaller respect to the errors in the pdf 's, but 
not insignificant for the rest. However, if one intends to take into considerations the 
impact of possible resummations on the pdf 's in a quantitative way, then the issues that 
we raise become crucial for obtaining a correct quantitative answer. In fact, it is also 
quite likely that, in the long run, with the large data flow from the LHC, these sources 
of errors that we investigate will become more significant in order to obtain more precise 
parameterizations of the pdf 's in the future. Our approach is part of an ongoing effort to 
develop a complete numerical program, Candia, that we hope will include not only an 
analysis of the evolution, with applications to DY and the Higgs sector through NNLO, 
with the inclusion of resummation effects, and so on. At this time Candia contains the 
DY and Higgs total cross sections beside, of course, the evolution. 

1.1 Comparison with the previous Uterature 

Regarding the constraints on our analysis, some comments are in order. In comparing 
the various results one possible source of disagreement lays in the treatment of the heavy 
flavours. Candia and PEGASUS treat the heavy flavours following the varying flavor num- 
ber scheme (VFN) as in [15j 0. MRST and Alekhin instead, follow different prescriptions, 
respectively described in [15] and [T7j. 

We have also found that other finer issues, in general not discussed in the previous 
literature, introduce systematic differences. For instance, MRST give a parameterized 
form for the input distributions at /Xq = 1 GeV^, while the lowest value of fij in their grid 
is 1.25 GeV^. In Alekhin's case, this author gives an analytic parameterization at 9 GeV^ 
without including the charm quark, even if the initial scale is above the charm threshold. 
The charm contribution is instead present in his evolved pdf's, available on a grid. The 
differences induced on the cross sections by the two methods are not negligible. A source 

^We will denote with the acronym VFN the varying flavour number scheme that follows the treatment 
of [H]. 
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of disagreement between Candia and Pegasus can be attributed to the fact that a given 
initial condition has to be fitted to a certain functional form in Mellin space, if one solves 
the equations using Mellin moments. These limitations are absent if one works directly 
in x-space [H]. 

1.2 Our approach 

As we have already mentioned, we base our analysis entirely on the implementation of an 
algorithm that solves the DGLAP equations directly in x-space and uses an ansatz based 
on various logarithmic expansions. These expansions have been shown to be related either 
to exact or to truncated solutions of the renormalization group equations (RGE's), char- 
acterized by coefficients which are determined recursively. Notice that these expansions 
are also typical of Mellin space [IH] , [S] • 

The structure of the recursion relations solved by this method is fixed by the choice of 
the original DGLAP equation and by the approximations performed on its right-hand-side 
(rhs), justified in a perturbative fashion. Beside the logarithmic expansions, also exact 
solutions are available for the nonsinglet sector up to NNLO, as we have shown, [H], which 
are useful in order to establish the convergence of the expansion toward the exact solution, 
whenever, of course, this is available. The method allows to bypass the appearance of 
commutators in the definition of the iterated solution from Mellin space, which is an 
unavoidable step in the singlet sector, as done in all the previous literature [T8] . [T9] . [20] . 
The logarithmic expansions, in both sectors, instead, take exactly the same simple form, 
being either scalar or matrix-valued. The relations for the unknown coefficient functions, 
introduced by the ansatz, and which are determined recursively, are also the same in both 
cases. 

In this paper we elaborate on the main features of these expansions by performing a 
thorough numerical analysis of the various truncated solutions introduced in our previous 
paper [H], discussing their behavior. In the singlet sector we also show the fast numerical 
convergence of the expansion and compare our results with those of the Les-Houches 
benchmarks which are based on a toy model of initial conditions. Comparisons are done 
both for a fixed and for a varying fiavor number. The anomalous dimensions involving 
the heavy fiavors have been implemented as in [15l [19]. 

We are going to see that variations induced by the choice of the solution induces 
variations on the cross section of the order of 1% or so at NNLO, and clearly affect also 
the NNLO ii"-factor for the total cross section. In our determination, the change in the 
value of the cross section from NLO to NNLO is around 4% on the Z peak, while the 
MRST and the Alekhin determinations are 2.6% and about 1.5% respectively. While these 
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variations appear to be more modest compared to the analogous ones at a lower order 
(which are of the order of 20% or so), they are nevertheless important for the discovery 
of extra neutral currents at large invariant mass of the lepton pair in DY, given the fast 
falling cross section at those large values. However, as we are going to show, the errors on 
the pdf 's induce percentile variations of the cross section as we move from NLO to NNLO 
of the order of 4% around the best-fit result, reducing the NNLO cross section compared 
to the NLO prediction and rendering these results compatible. 

2 Initial sets 

Our comparisons are performed using an implementation of the theory in a code written 
by us and called Candia that will be documented elsewhere. The other implementation 
that we can directly compare to is the code written by A. Vogt, Pegasus [I9], which 
implements the Mellin-transform method. This is the only NNLO code which is of public 
domain at this time. Pegasus can be run in different modes and allows to select numerical 
solutions of a given accuracy. Our evolution is also compared to the MRST evolution. We 
clarify in our results if we have used in our implementations of the hard scatterings the 
MRST input evolved by MRST or our evolution of the same input. For the evolution of 
the MRST parametric input with Candia we have worked in the same VFN scheme, but 
we have used a slightly different prescription, and the comparisons are performed either 
starting from their values on the grids or from the parametric input provided by the same 
authors. 

As we will show below, the grid and the parametric inputs (the first at /Xq = 1.25 
GeV^, the second at /ig = 1 GeV^) give results which differ at the percent level. 

The closeness between our predictions, obtained using the asymptotic solution, whose 
nature we clarify below, and their distributions is also quite evident, but the variations 
are such to generate differences at the percent level in the cross sections. Candia gener- 
ates numerical outputs of the exact solutions up to NNLO for the nonsinglet sector and 
truncated solutions with arbitrary order of accuracy in the singlet sector. In our case the 
term "arbitrary order of accuracy", referred to a solution, means that we use logarith- 
mic expansions in log(as) multiplied by coefficients of a certain power of as rather than 
log(x(«s)), where x{(^s) is a typical NLO function of the running coupling or some other 
non-trivial (composed) function generated at higher orders. We remark that the simple 
\og{as) expansion converges after very few iterates (7-8), with a precision of 4 to 5 signif- 
icant digits. Clearly, if one is searching for exact solutions, the iterates converge rather 
slowly to give the "exact" numerical solution since the leading order solution is not fac- 
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tored out in the ansatz, as in the case of the [/-ansatz of [T8l fT9]. as we will clarify below. 
The differences in the cross sections are very small (0.1 %), if an asymptotic truncated 
solution (with k' = 7 or 8) replaces the "exact" solution, or brute force solution. But the 
theoretical indetermination remains: at NNLO even the second truncated solution is a 
solution and the differences on the observables, as we are going to show, become more 
substantial than the fraction of a percent obtained using the asymptotic solution. 

As we are going to show next, the results produced by our implementation are in 
excellent agreement up to NNLO with those obtained with Pegasus and the Les Houches 
benchmarks (see refs. [21], [22] )• Regarding the computation of the errors on the pdf's, 
these are not available for all the most popular sets and through all orders. For instance, 
the Alekhin and MRST fits are presented up to NNLO, but only one of them, the set [9] 
presents errors through NNLO. 

Coming to describe the region that we have studied, our numerical investigation covers 
both the resonance region around the peak of the Z - this being useful in order to assess the 
impact of the corrections at the various orders - and the remaining regions of faster fall-off. 
We remark that these studies can also be extended to the search for extra Z' in extensions 
of the Standard Model, and as such provide a clear indication of the role played by these 
corrections in a precise determination of the QCD/electroweak background, useful for 
potential discoveries of additional neutral currents at the LHC using this process. 

2.1 A classification of the possible solutions 

We start our study by briefly reviewing the nature of the NNLO solutions and the level 
of accuracy which is intrinsic to any solution. 

As we have already mentioned, there are essentially three types of solutions that one 
can extract from the evolution equations. We have decided to classify them as follows: 1) 
the accurate solutions with few iterates, 2) the iterated solutions with a large number of 
iterates, eventually combined with exact analytical solutions of the nonsinglet sector and 
3) the brute force solutions. 

Solutions of type-1 are built using a very simple ansatz which can be showed to be 
accurate through NNLO. The ansatz contains all the terms of the form log"'(Q;), log"(a) 
and al log" (a). Solutions of type-2 include terms of arbitrary higher order in the logarith- 
mic expansion. They are characterized by one or two indices, the first (k') which defines 
the accuracy and the index k which defines the order of the truncation of the right-hand- 
side (rhs) of the evolution equation. The two indices be taken as label of a given 
truncated solution. They are built starting from one of the two forms of the evolution 
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equations (form-1 and form-2), as we will discuss below. The index k appears in form-2, 
when the DGLAP equations are written directly in terms of the logarithmic derivative of 
the running coupling log(Q;s). 

Solutions of type-3 are not accurate since they are affected by contributions of all 
orders and can be obtained using brute force methods, by discretization of the evolution 
equations. They exceed the level of accuracy typical of a perturbative expansions where 
the kernels are only known up to NNLO (a^). 

Truncated solutions, instead, are obtained, as we have just mentioned, by truncating 
the rhs of the evolutions equations at a given order and then searching for solutions of 
these equations with a certain accuracy. Also in this case the truncated equations can be 
solved exactly, giving solutions which exceed the level of accuracy of the expansion. This 
happens in the nonsinglet sector. If we use the DGLAP written in terms of a logarithmic 
derivative of Q rather than of the coupling (form-1), then the NNLO exact nonsinglet 
solution is available. Similar solutions are also available for the form-2 of the equations in 
the same sector for k, < 2. In this second case these solutions do not have, however, a well 
defined meaning, since they are not accurate nor converge, in the limit k' oo, to the 
exact solution of the exact equation (as for the ansatze written for form-1, for instance). 
Therefore, also in this case an accurate solution is obtained by an additional expansion 
in the strong coupling and retaining only terms up to a certain order, which is the order 
of the selected accuracy. 

2.2 The two forms of the evolution 

We proceed by illustrating the two forms that the equations can take. 

We recall that the perturbative expansion of the DGLAP splitting functions and of 
the /9-function take the generic form (to the m-th order) 



with /?fc being the corresponding coefficients of the f3 function which have been summa- 
rized in [ll]. Leading, next to leading and NNLO correspond to the cases m = 0, 1, 2 
respectively. 

The equation can be written either as 



P^'"^°(a;,Q2) = ^('^£g!)') 



P«(x,Q2) 





df{x,asm) 
91og(g2) 



P 



,N'^LO 



{x,Q')(E)f{x,a,{Q')) 



(2) 
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(form-1) or, equivalently, as 



d fix, as) P^"^^Q(a:,g2) ^,, 

- pNm,o ^f(^^^s) (3) 

(form-2). While the two forms are equivalent, form-2 needs an expansion of the 1/(3 
factor. This generates on the rhs of the expanded equation an infinite set of truncated 
equations characterized by a parameter of accuracy (/t). This parameter has to be sent to 
infinity in order for form-2 to be equivalent to form-1. When we search for solutions of the 
DGLAP in the form-1 and we need to compare with the form-2, the recursion relations 
for the solution ansatz start to differ after the order k. Notice that once that we have 
introduced an expansion of the rhs, such as in form-2, we may search either for truncated 
solutions of this truncated equation or for exact solutions. These are options that increase 
the type of possible solutions, all of them of different theoretical accuracy. To illustrate 
this point we start from the form-2 of the equations and choose m — 2 (NNLO), obtaining 
the truncated equation 

df((x Xi 1 

gj' ' = — [Ro{x) + a^Riix) + alR2{x) + • • • + ® /{a^, x). (4) 

The explicit form of the operators can be easily worked out, at any perturbative order. 
For instance at NNLO the first few terms are given by 

Ro = 

R^ = - J-P(i) - -\p(o) 
7r/3o (47r) 

R^ — — ^r-^^^^ ~ 7j~\oRo ~ / \ Ri 

-^4 = — J-. — ^Ri + — r^R2 — J-. — ^R2 
(47r)3 (47r)2 (47r)2 

bl bl 6162 



(5) 

where 6j = Pi/ Po, valid in the nonsinglet case. A similar expansion holds also for the 
singlet sector, although in this case the recursion relations involve some commutators of 
the matrix-valued kernels. The unknown operators that define the ansatz need to be 
identified by solving the related recursion relations. 
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In Mellin space, the ansatz that solves the Eq. (jll), is chosen to be of the form [TSf fT9] 



U{N, as)fLo{N, ao)U-\N, «o) 



+00 



k'=1 



fLo{N,as,ao) 



+00 

k'=1 



-1 



(6) 



that we call, for convenience, the U -ansatz, and inserting this expression in Eq. (jl]) it 
generates a chain of recursion relations which allow us to determine the matrices Uk>{N) 
in terms of the operators -R^- 

We remark that even if we take a fixed value of k in the truncated equation, the 
running index k' in Eq. is still free. The case k' > k with — > 00 § allows to find the 
exact solution of the K-truncated equation (jlj). 

A third option corresponds to the choice k' < k. This gives an approximate solution 
of the K-truncated NNLO equation accurate at the order af. Notice that if we start from 
the first form of the evolution (form-1) and use a recursive ansatz to solve this equation 
(either in moment space or in x-space) this solution has to agree with the solutions of the 
truncated equation considered above, once we perform an expansion of that solution in 
as and ao, as we have showed in 



As an example, the most accurate NNLO solution is generated by the choice n' = k = 
m = 2. In this case we can write the a^-truncated solution of the truncated equation 
with K = 2 as follows 



f{N,as 



"0 



.J-p(O) 



[l + (a,-ao)f/i(iV) + «'f/2(A^) 
-asa,Ul{N) + al {Ul{N) - U^iN))] f{N, ao) 



(7) 



At this retained accuracy (m = 2) of the evolution integral, the truncated solution of 
the corresponding (truncated) DGLAP equation can be easily found, in moment space. 



as 



fiN,as) = f{N,ao) 



a. 



p(0) 



1 + (a. 



ao) 



p(i) p{o)p^ 



+ai 



+"0 



par p(2) 



2772/32 47r2/5o 



pWpWp^ PWP^ p(0)2^2 p(0)^2 p(0)^^ 



2vr2/32 



^vr2/32 



167r2/33 



p(i)^ 



+ 



p(2) 



p(o)p(i)^^ p(i)/3i p^^rpl p(o)(3f 



27r2/52 47r2/5o 27i^(3^ 8n^(3? 







167r2/32 



^ This corresponds to the option IMODEV = 2 in Pegasus while the with 
corresponds to the option IMODEV = 1. 
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These are solutions of type-1. They coincide with the first few terms of the exact NNLO 
solution of the DGLAP equation, obtained by a double expansion in the couplings and 
retaining only the 0(q;^) terms, as can be explicitly checked in the nonsinglet sector for 
the equation given in form-1. Therefore the solution is organized effectively as a double 
expansion in a.^ and ao- This approach remains valid also in the singlet case, when the 
equations assume a matrix form, though an exact solution, in the form of an ansatz, 
similar to that of the nonsinglet sector (Eq. [HI] below), is not available in this case. 




2.3 The logarithmic expansions for the form-1 of the evolution 

Our previous analysis has involved form-2 of the equations and we have presented an 
ansatz that solves this equation. We intend now to show how to construct an ansatz 
directly starting from form-1. 

The advantage of solving the equations directly in form-1 is that one has a single 
ansatz for the entire equation and the accuracy is just determined by the order of the 
chosen ansatz, differently from form-2. If we are interested in an accurate solution of 
order k', for instance, we use the ansatz 

oo 
n=0 

(9) 

which can be correctly defined to be a truncated solution of order k' of the DGLAP in 
form-1. As we are going to show next, we will monitor the numerical behavior of this 
expansion and its convergence. Sending the index in the logarithmic expansion of ([9]) 
to infinity, then the ansatz that accompanies this choice becomes 

(10) 

and converges, in principle, to the exact solution of the equation given in form-1. In 
practice, however, this convergence is hampered by the factorial suppression. For this 
reason is it convenient to use the term "asymptotic solutions" rather than "exact solution" 
for these iterates of larger index. 



X] 



_n \ 1 n 



In 



c^sm 
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2.3.1 Exact Solutions in the nonsinglet case at NNLO 



The search for exact NNLO solutions in the nonsinglet sector proceeds similarly. This 
has been analyzed in [ll]. We define the following functions 

>C = log-, (11) 

°^167r2/5o + 47rao/3i + a2/52' ^ ' 

.(iV) ^ (14) 

W = — (15) 

c(iV) = Mlp(o)(Ar)-8P«(iV) + ^p(2)(^)^ ^^Q) 

Po P2 

where for n/ = 6 the solution has a branch point since 4/3o/?2 — /^i < 0. If we increase as 
we step up in the factorization scale, for n/ = 6 Q is replaced by its analytic continuation 

Q = — =arctanhY (17) 
VP! - 4/3o/32 

where 



27r(a,-ao)\/W2-/5? 
X = 

27r(87r/3o + (a^ + ao)/5i) + asao/?2 

1 , /I + X 
arctanhx = - log 



2 Vl-X. 

Clearly all the (nontrivial) dependence on the coupling constants is contained in the 
3 functions C,J^ and Q. The general solution can be written in terms of A'^,B'^,C'^, 
coefficients that will be calculated by a chain of recursion relations [11] giving 

fM') = (e^':") (E^-W'") (E^S') nx,Ql) 

\n=0 ■ / (g, \m=Q ■ / \p=0 ^' / 

OO S t j~\g / \ 

t^t^£^^'(t-n)!(.-t)! 



and where 



DtJx) = <(x) ® P^„(x) ® C;_,(x) ® fix, Ql). (20) 
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Solving the chain of recursion relations, the above solution in x-space can be simply 
written as 



f{x,as{Q^)) = exp 



exp \ log 



exp 



/ 167r^/3o + 47ra,/3i + ag/^a 
V 167r2/3o + Anaopi + ag/^s 
1 27r(a 



p(0)| 



x) 



4p{2), 



X) 



f32 



arctan 



2^: 



P2 } ) ' 



(21) 



where D, 



0,0 



X 



/(x, Qo)- The possibility of finding an exact solution has, of course, 
phenomenological implications, since the analytic solution performs a resummation of the 
log(Q;s) which are generated to all orders by the various truncations and by the corre- 
sponding logarithmic expansions. These are incorporated into the functions Ai, Q (x). 

We will get back to this point later. 



2.4 The Singlet Case 

Before we address the topic of the resummation/re-organization of the logarithmic struc- 
ture of the solution due to the choice of the different expansions, we move to analyze the 
extension of our previous reasonings to the singlet case. One can start from form-1 or 
from form-2, obtaining solutions of overall different accuracies. In the singlet case, if we 
start from form-2, then one can consider a truncation of this equation, for instance to 
second order, that can be written as 



where 



df{N,a,) _ 1 
da, a. 



Rq + agRi + Oi R2 



(22) 



Rq — 

Ri^ 
4 = 



A 



.p{0) 



^ R^ ^ i?o/52 



(23) 
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whose (exact) solution in Mellin space is expected to be of the form (the [/-ansatz) [18] 



f{N,as)= I + asUi{N) + alij2{N) L{as,ao,N) 1 + «ot>i(iV) + «of>2(iV) f{N,a^) 



(24) 



where 



-R2 - RiUi + 2U2. 



(25) 



Using two projectors on the subspaces of the corresponding leading order (singlet) eigen- 
values, denoted by {e±) (see [E]), one can remove the commutators, obtaining 



where 



U2^ 



U2 = Ut^ + Ut + f/2-+ + f/2-- , 



+ Rt+ ^1 



(26) 



R]^ R^ + i?2 
1 



r_ — r_|_ — 1 



r_|_ — r_ — 2 
1 

r_ — r I — 2 



r_|_ — r_ — 1 
Rl-R{- - R+- + 



r+ — r_ — 1 
r_ — r 1 — 1 



(27) 



and the formal solution from Mellin space can be simplified to 



L + ttsUiL — aoLUi 

- asaoUiWi + all (u^ - U2) /(iV,ao) (28) 



f{N,as 



where the accuracy is kept through 0{al). 

If we don't want to truncate the equation, then we work with form-1. We start 
constructing solutions of this equation using the logarithmic expansions introduced in [H] 
using few iterates. As we have mentioned, in this case there will be just one parameter 
appearing in the expansion, related to the desired accuracy, i.e. the terms retained in 
the ansatz, and by increasing the accuracy one expects the result to converge toward the 
exact solution. 
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The first truncated logarithmic ansatz that is expected to reproduce (l28l) includes also 
an infinite set of new coefficients C„, similar to the nonsinglet NNLO case 



f{N, as) = J2— [An + a,Bn + a^Cn] . (29) 



n! 

n=0 



The ansatz can be generated to an arbitrarily high order. If this order is k, we introduce 
the K-truncated logarithmic ansatz 



In 



n=0 

in the NNLO DGLAP matrix equation and neglect the 0(0;^+^) terms. We obtain the 
following recursion relations which, in the NLO DGLAP case are 

0°+i(x) = --^ [P(°)(a:)®0°(x)], 
Po 



(30) 



-Ol-[{x) - kO:{x) -{k- l)£^or^(x) , (31) 



while in the NNLO case become 

0°+i(x) = -|- [pW(x)®0°(x 
Po 

Ol^,{x) = --^ [pW(x) ® Olix)] - [P(i)(x) ® 0°(x 
Po ttPo 



47r/3o 



0°+i(a:)-0^(a;), 



0:;+i(x) = -|- [P(°)(x) ® 0:;(x)] - ^ [P(i)(a:) ® Ori(x)] 
Po Trpo 



L_ rp(2), 



x) ® or^fx 







-nO^ix) -in- l)^0:-\x) - (n - 2)-f^OT'{x) . (32) 

These relations hold both in the nonsinglet and singlet cases and they can be solved 
in x-space and N-space in terms of the initial conditions f{N, a^) = Oo(A^). Since in the 
singlet sector the recursion relations are in matrix form, we can solve them by the use of 
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the projectors e+ and e_. A straightforward way to solve the matrix relations is first to 
solve the relation for 0° (A^) in terms of e+ , e_ and r" , r" as follows 

0° (AT) = e+r!^0°(Ar) + e_r!!0°(Ar) (33) 

and use this result to solve the other relations. The 0™(A^) operators can be decomposed 
in an orthonormal basis {61,62} as 

0^(N) = eiO«'"^(Ar) + e20i')'"^(A^) = 0+'"^(A^) + Q-'^'iN) . (34) 

Then, using the properties of the projectors we can write 

(e+ + e_) O-(Ar) = e+ (eiOW^'"(Ar)) + e+ (620(f)'-(Ar)) 

+e_ (6iOW'-(Ar)) + e_ (620(2)'- (A^)) , (35) 

for the relations with m > 0, and setting 

0++'"^(Ar) = e+ (6iO«'"^(A^)) 
0+-'-(Ar) = e+ (620(f)'-(Ar)) 
0;+'-(Ar) = e_ (6iOW'-(Ar)) 

0;-'-(Ar) = e_ (e20(f)'-(Ar)) (36) 

we can derive some recursion relations. For example, in the NLO case, which corresponds 
to the case m = 0, 1, we have two recursion relations and having solved the 0°(A^) as 
illustrated above, the m = 1 relation can be decomposed into four recursion relations as 
follows 

l)0--'\N). (37) 



0::i\N) = R++r!^0°(Ar) + (r+ - 
Oti\N) = R+-r!^0°(Ar) + (r_ - 
0-ti\N) = Rr+r^O°(Ar) + (r+ - 
Otli\N) = RrrlO',{N) + (r_ - 



This pattern can be extended to NNLO, in fact we have three sets of recursion relations 
corresponding to the cases m = 0, 1, 2. Once we have solved all the relations corresponding 
to the cases m = 0, 1 we can proceed to solve the following relations 



(AT) 



R++ + 



01 



/3o(47r) 



/5o(47r) 



R++ rlO',{N) + [R++0:+'^(Ar) + R^ 0;+'^(Ar)] 
[O^-'^'iN) + 0::i\N)] + -^^r+0:^'\N) + (r+ - 2) O^^'^ , 
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Po(4vr) 



-"0°(iV) + \^X-0--\N) + R++0+-'i(iV)] 



,2 

n 1 



Po(4vr) 



-"o°(iv) + [Rr+o++'^(iv) + Rr-o;+'^(iv 



Pi 



Po(47r) 



^!^0°(iV) + [Rr-0;-'^(iV) + Rr+0+-'i(iV)] 



Po(47r) Po(47r) 



(3^ 



which can be implemented in a computer program, with a standard numerical inversion 
of the Mellin transform, being equivalent to (l3Ti) and (1321) . The x-space approach, as we 
are going to show, matches the numerical Mellin method with very high accuracy, since 
the asymptotic truncated solutions give the same answer. In the nonsinglet sector the 
exact solutions built by iterations as logarithms of composite functions of are new and 
not present in the previous literature. These have been used in this sector to generate the 
corresponding exact solutions. 



2.5 Relating the [/-ansatz to the logarithmic expansion 

It is important to compare the two expansions which are identical globally (that is to all 
orders) but that organize, at a certain fixed perturbative order, the logarithmic corrections 
in different ways. This can be easily shown in the nonsinglet sector, where the two 
expansions can be more easily mapped into one another. Let's see how this happens. 

The double Taylor-expansion of the solution of the Eq. (12^ for {as,ao) around (0,0) 
up to order 4, for example, has the following structure 



f{x,as,ao) = { — 
ao 



1 + agtt^i^ + aoa^^ + 



2 (2) , (2) , 2 (2) , 
a^al + asaoa2 + aoO-s + 

3 (3) , 2 (3) , 2 (3) , 3 (3) , 2 2 (4) 

4 (4) , 3 (4) , 2 2 (4) , 3 (4) , 4 (4) , ,44 (8)' 



/o , 
(39) 
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as we can see, the double expansion gives terms of higher order of the type a^aQ. On the 
other end, for instance, the logarithmic expansion accurate to is given by 



f{x,as,ao) = \An{x) + asBn{x) + a^C„(x) + alDn{x) + a^E„(x)l ^log" ( — ) 



(40) 



that gives recursion relations for the coefficients An-, ■ ■ ■ ■, En which are solved and expo- 
nentiated, as we have shown in |T1]. Once those coefficients have been determined, we 
substitute them into Eq. (I40l) and rewrite / as 

fa \ ~^^^°^ 

2 (2) , (2) , 2 (2) , 

a^Ci + a^aoCj + a^c-^ + 

3 (3) , 2 (3) , 2 (3) , 3 (3) , 

ajc^^^ + asalc^2^ + ala^c^^^ + aga^c^^^ + a^c^^^ 



Ao. (41) 



From the direct calculation of the coefficients and c^-"'^ in the two Eqs. (HTl) and ( l39ll . 
we observe that they coincide only for those terms which are of the same order in a^ao, 
but in general, the two expansions organize the corrections in different ways. For instance, 
in order to generate the terms of the type a^aQ, we should take the index k = 8 in ( !30l) . 
In this case we will reproduce all the coefficients up to a^a^, but we will also introduce 
terms of order and a^ which were not present in the double expansion of (12^ arrested 
at order a^. This is due to the fact that the Taylor expansion of in (a^, Oq) is a double 
expansion while the result of the logarithmic expansion corresponds to a single expansion 
in as and the remaining power of ao are introduced during the exponentiation procedure 
[H]. As we have mentioned, to establish the equivalence between the two approaches 
Eqs. fl^ and flHUl) one needs to expand the leading order solution which appears as first 
factor in fl2^ . extracting all the logarithms of a^. The structure of the f/-ansatz is such 
that in it the leading order solution is automatically factored out, while in the logarithmic 
expansions of type (129|) and, in general, ( l30i) . one needs to exponentiate the solution of the 
recursion relations to achieve the same result. Numerically this can't be done, but the two 
ansatze, interpreted perturbatively both as ways to collect the logarithms of the solution 
of the evolution equations, become the same expansion as the order of the truncation 
grows. 
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3 Resummation and the exact solution 



It is interesting to compare the logarithmic corrections generated by the truncated solu- 
tions with the exact nonsinglet solutions obtained at the various perturbative orders. As 
we have already mentioned, the analytic solution resums the partial contributions coming 
from the truncates of various order introduced by the various ansatze in x-space or in 
moment space. To illustrate this point, let's start the analysis from the NLO nonsinglet 
case and then we will generalize the results to the NNLO case. 

Solving NLO DGLAP nonsinglet equation in Mellin space 



dots PNLOyOis) 

we obtain an exact solution which can be written as follows 



(42) 



/(«.,iV) 



exp{-lp(o)^iV)log(^)} 



X exp 

in Mellin space, and as 



,Po Pi 



log 



47r/3o + UsPi 
47r/3o + QioA 



/(ao,iV) (43) 



fM = exp{-lp(o)^,),,^^^^| 



X exp 



|-p(°)(x)--ipW(x) 

,Po Pi 



log 



A'k(3o + asf3i 



47r/3o + aoA 



/(ao,x) (44) 



m x-space. 



Expanding in terms of log this solution we obtain 
f{as,x) = exp|--^p(°)(x)log 



exp 



|-P(°)(x)--ip«(a;) 
Po Pi 



47r/3o + q;o/3i 



a 



log - + 



/(ao), 



(45) 



with an analogous expression in moment space. The notations can be simplified by defin- 
ing 

p(0) = _Ap(o) 

p(i) ^ A p(o) _ £ p{i) 

/3o Pi 



5'i(ao) 
t = log 



"0, 



(46) 
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and in x-space we can rewrite the solution in terms of t-iterates in the form 

exp|p(°)t| exp|p«t(7i(ao)| exp jp^^) ^2(«o)| ■■• 
exp |p(^)rc/„(ao)| /(ao,a;) 



where 



^2(«o) = ^ (^i(«o) -^?(tto)) 
^3(ao) = Q^i(ao) - ^^?(ao) + ^^?("o)^ 

£/4(ao) = (^^i("o) - ^c/i(ao) + ^^i(ao) - ^gti^o)^ 
Finally, in the nonsinglet case we can re-arrange our solution in the form 



(47) 



(48) 



fLo{<ys,x) = exp \P^°h\ f{ao,x) 

/(a„x) =exp|pWtc/i(ao)| exp jp^ ^2(^0) | ■ ■ -exp jp^^) c/„(ao)| /Lo(as,a;). 



(49) 



It is interesting to note that the function gi{ao) is, in a sense, universal since it contains 
all the information about the initial conditions. A quick comparison between (jH]) and its 
expanded version ( l49l) shows the features of the implicit resummation involved in mov- 
ing from the second equation to the first. We will point out, in the numerical analysis 
presented below, that only a resummation can bring a logarithmic ansatz expressed in 
terms of log(as) (either in Mellin space or in x-space) to reproduce numerically the ex- 
act solution. This is easy to show in the nonsinglet case, where both equations can be 
implemented as numerical iterations. 

In a similar way we can proceed to re-arrange the exact solution in the nonsinglet 
sector at NNLO. This can be rewritten as 



f{x,as 



exp < log 



pW(x) _ 4P(')(x) 
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exp 



^pW(x)-8PW(a;) + ^p(2)( 



arctan 



27r(87r/3o + {oig + q;o)/3i) + q;sQ;o/32 



X] 



(50) 



Expanding in terms of the logs, it is useful to define the following expressions 

G{ao) = 



(51) 



Then we get 



/(a;,a.) ^ exp {^a,(ao)Pm ^ 



a2(ao)^A ^ + h ^„^^^-^ an(«o)^A^ 



exp 



^(cvo) 



&2(ao)P?^ + --- + 



^^(cvo) 



&n(ao)^B^ fLo{x,ao) 



(52) 



where ^(ao) and the functions ai(Q;o), . . . , &i(q;o) • • •, are polynomial functions dependent 
on ckq. We omit to give their explicit expressions since they are not relevant for our 
discussion. With these definitions, the solution written in terms of simple logarithms of 
the coupling is summarized in x-space by the formal expression 



f{x, as) ~ exp 



exp 



ai{ao)P^^\ exp ^ 



a2{ao)PX 



(2) 



-bi{ao)PB^ } exp 



G{aQ 



G\ao) 



b2iao)P^B^ 



I exp 



bn{ao)PB^ \ fLo{x,ao). 



G-{ao) 



■ exp 



G-{ao) 



(InMPA 



(5 



The relations between exact solutions and logarithmic expansions simplify consider- 
ably when one starts from the form-2 of the evolution equations. In fact, proceeding with 
the 1st truncated equation {k — 1) this takes the form 



dficx X) 1 
•'^ ^' ^ = —[R^ + asRi]^f{as,x), 



(54) 
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where we have set 

i?o = -|-P(°) i?i = -P«^ + pW-^. (55) 
In this specific case the exact solution is given by 

f{a,, x) = exp {{as - ao) ^i}^^ exp {tRo}^ /(ao, x) (56) 
and using the relation 

f{as,x) = exp{tRo}^exp{aotRi}^exp l^ao^—Ri'^ ■■■f{ao,x) (57) 

followed by a further expansion of the exponentials, the expression above can be re- 
organized in the form 

/(a„x) = exp{tRo} ^ ll + Riaot + f ( Ri^ + Ri ^ Ri^") + ■ ■ ■] (g) f{ao, x). (58) 



If we want to preserve a certain accuracy in our solutions, it is sufficient to do a Taylor 
expansion of fl3U|l . For example, at NLO, the truncated solutions of order of the 
truncated equation is 

f{as, x) = [l + (as - ao) Ri] (g) fLoi^s, x), (59) 

which takes the form originally given in [20] • Expanding this expression in terms of 
log(as/ao) = t we obtain the traditional form of the solution 



f{as,x) = fioiots^x) + Ri 



1 2 



fLo{as,x). (60) 



Using this simple approach we can proceed to the determination of finite accuracy 0(a;^) 
solutions in the nonsinglet sector. 

Increasing the value of k, we can write the K-th truncated NLO or NNLO equation as 

o f ( (y X) 1 

^J' ' = — [Ro + asRi + alR2 + ... + a^R^] ® /(«„ x), (61) 

where all the coefficients Ro{x), Ri{x), . . . , Rf^{x) are expressed in terms of the P^^^ and 
P(i) kernels in the NLO case, and in terms of P^^\ P^^\ P^^'' in the NNLO case. In both 
cases the solution can be expanded in terms of t-logs as 

f{as, x) = exp {tRo} exp [t (ao-Ric} + alR2c\ + . . . + aoP^c;^) } ® 
exp [t^ {aoRicl + alRicl + ... + Oo^kC^) } ® • • • 
®exp {r (ao^iCi + ao^2C^ + . • . + a^P«0 } ® /(ao,a;), (62) 
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being the coefficients real numbers. After a further expansion one can cast the result 
in the form 

f{as, x) = {l + t [aoRicl + alR2c\ + . . . + a^R^^c^^ 

[aoRicl + alR2cl + ... + a^R^cl) ® 
(ao-RiCi + alR2cl + . . . + ao-R^c^) ® ■ ■ ■} ® fhoiao, x), 

(63) 

having factored out the leading order solution. 

One of the points that should be briefly taken into considerations concern the defini- 
tion of the asymptotic solution. An asymptotic solution, in our terminology, identifies a 
solution which is the closest possible to the exact (brute force) solution. This means that 
while in the nonsinglet, for this solution, we will be using our exact ansatz, for the singlet 
we will let the number of iterates grow until the logarithmic series stabilizes. However, the 
absence of exact solutions in the singlet case shows that we will be surely differing from 
the brute force solution by some finite amount. Being Candia, or Pegasus based on 
analytical approaches rather than on discretizations, we are not able to compare with the 
exact solution and estimate the difference between our asymptotic solution and the exact 
one. We will quantify these difference rather accurately taking the Drell-Yan cross section 
as an example, but before coming to a numerical analysis we discuss the implementation 
of the renormalization scale dependence in our formalism. 



3.1 The treatment of the renormaUzation scale dependence and 
the implementation 

The scale dependence of the pdf 's can be obtained by solving the modified equations 
d 

fi{x, /i^) = Pij{x, /i^) ® fj{x, fx%) , (64) 



d In fij^ 

where /Zi? is now a generic factorization scale. The explicit expression of these modified 
kernels are given below [19]. This can be obtained by re-expressing the coupling constant, 
function of the factorization scale fip, in terms of /i/j using the RGE for the running cou- 
pling at the corresponding order. Concerning the actual relation between the couplings at 
the two scales, this can be obtained by solving numerically the corresponding RGE for the 
running coupling at NLO and NNLO. We have also monitored the approximate solutions 
obtained by the usual well-known asymptotic expansions in terms of L = ln(/x|,///|j). In 
the NLO case an implicit solution which allows to connect and /i^ is available 

' - ' +/?olnf4V&iln|^44^j^4##| (65) 



asif^p) a^(/^|) Xf^R/ las(/Lt|) [1 + 6ias(/i|,)] 
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where as{fj,'^) = cts (/x^)/(47r), which can be solved as 



2/ 2^PoL , 



(66) 



where the dependence is contained in the factor L, and we have used a /3-function 
expanded up to NLO, involving Pq and At NNLO implicit solutions such as are 
not available but one can derive the analogous of ( |66l) . Both options, the exact and the 
asymptotic are present in Candia. The differences between the two determinations are 
quite small (see Tab 1201) . 

We have imposed logarithmic expansions on the equations with the kernels written in 
the form given below and derived recursion relations for these expressions. These reduce 
to the recursion relations discussed in the previous sections with the actual redefinitions 



«s(/iR) = as(/i|) 



^2f,2\f^oL «s(A''f)/ n2T2 n t\ 



(47r)2 

^p(O) 

2 



p!--Hx) 



(0), 



x] 



(67) 



introduced into the equation expressed in form-1. 

Concerning the implementation of the algorithm in Candia, we briefly illustrate the 
implementation of the flavor reconstruction. We define 

q,±q,, gW = ^gf), (6^ 
1=1 



(±) 



then the general structure of the nonsinglet splitting functions is given by 



P = P-- 



5,kPl + P^ 



qqi 



(69) 



Pq^qk - P^^<lk - ^ikPqq + Pqq- (70) 

This leads to three independently evolving types of nonsinglet distributions: the evo- 
lution of the flavor asymmetries 

A%. = (71) 
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whose evolution is governed by 



and the sum of the valence distributions of all flavors q^^^ which evolves with 

PnS = Pqq - Pqq + ^/ {Pqq " P^^ = PnS + PnS' (73) 

Notice that the quark-quark splitting function Pqq can be expressed as 

Pqq = Pns + nf [P^q + P'qq) ^ P^s + Pps. (74) 

It is important to observe that the nonsinglet contribution is the most relevant one in 
Eq. (1741) at large x, where the pure singlet term Pp^ = P^^ + P^^ is very small. At small x, 
on the other hand, the latter contribution takes over, as xPps does not vanish for x — 0, 
unlike xP^g. The gluon-quark and quark-gluon entries are given by 

Pqg = rifPq^g, (75) 

Pgq = Pm. (76) 

in terms of the flavor- independent splitting functions Pq^g = Pq-g and Pgq- = Pgg.. With 
the exception of the first order part of Pgg, neither of the quantities xPgg, xPgg and xPgg 
vanish for x —>■ 0. 

In the expansion in powers of of the evolution equations, the flavor- diagonal (va- 
lence) quantity P^ is of order as, while Pj^ and the flavor- independent (sea) contributions 
Pqg and Pqq are of order a^. A non- vanishing difference Pgg — Pgg is present at order a^. 

The next step is to choose a proper basis of nonsinglet distributions that allows us 
to reconstruct, through linear combinations, the distribution of each parton. The singlet 
evolution gives us 2 distributions, g and q^^\ so we need to evolve 2nf — 1 independent 
nonsinglet distributions. We choose 

1. q^^\ evolving with -P/J^g; 

2- Qn-sm = Qi'^ - Qi'^ (for "2 <i < rif), evolving with P^^, 

3- Qn-sm = Qi^^ - Qi^^ (for 2 < i < %), evolving with P^^, 

and use simple relations such as 



■I \ k=l,k^i / 



(77) 
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Figure 1: Plot of the percentage differences between the asymptotic Drell-Yan cross 
section and those obtained using expansions of the pdf 's of a fixed order k, shown as a 
function of k for the NLO and NNLO cases. We have used the MRST parametric input 
with fiQ = 1 GeV and Q = Mz- The evolution is based on Candia. 



to perform the reconstructions of the various flavors. Choosing i = 1 in ( 1771) . we compute 
q[ ■* from the evolved nonsinglets of type 1 and 2 and q[~^^ from the evolved singlet g^^-* 
and nonsinglet of type 3. Then from the nonsinglets 2 and 3 we compute respectively ■* 
and ql^^ for each i such that 2 < i < Uf, and finally gj and g,. 



Moving from NNLO to NLO things simplify, as we have 



5,(1) 

qq 



Pq/^'^ ■ This implies 



>v,(i) 



Pj^s^\ i.e. the nonsinglets q^ ^ and gjvsiit evolve with the 



(see Eq. ([7H])) that P^cj 
same kernel, and the same does each linear combination thereof, in particular q^~^ for 
each flavor i. The basis of the 2n/ — 1 nonsinglet distributions that we choose to evolve 
at NLO is 



1. g,[ ^ (for each i < Uf), evolving with P^ '^^^ 



NS ' 



2. q^^gii = q[^^ — g-'^'* (for each i such that 2 <i <nf), evolving with Pps^\ 



and the same we do at LO, where we have in addition Pj^g 



(0) 



P^'P, being Pgg 



v,(o) 
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4 The cross section and the parton luminosities 



Our NNLO analysis of the total cross section for lepton pair production combines the hard 
scatterings of [12], implemented by us in a program called CANDlA^iy, which combines the 
hard scatterings with the evolution performed by Candia. We will present in a section 
below some results obtained by interfacing Vrap and Candia that allow to extend some 
of the predictions of [13] with the inclusion of the factorization/renormalization scale 
dependence not only in the hard scatterings but also in the evolution. Here our main 
analysis is instead focused on the cross section for the mass distribution da/dQ"^. 

Lepton pair production at low Q via the Drell-Yan is sensitive to the pdf's at small- 
X values while in the high mass region, above the peak, is essential for the search of 
additional neutral currents. The general structure of the factorization formula for the 
color averaged inclusive cross section for lepton pair production is given by [12] 



dQ^ 

where ay is the point-like cross section in the case of the 7, Z and the interference 7-Z. 
S is the center of mass energy of the incoming hadrons and is the invariant mass of 
the di-lepton pair, respectively. We have used the relations 



oz{Q 



"em 



2^ 



3g^ Nc 

AMzNc sin^ Ow cos^ Ow {Q^ - M|)2 + M|r| 
TraL (l-4sin^M (Q' - M|) 

6 sin^ 9w cos2 Ow NcQ\Q^ - M^f + M|r| ' 



(79) 



where T z^u = 0.0839136 GeV, Tz = 2.4952 GeV, sin^ 9w = 0.23143 and a^UMz) = 
1/128. These choices, performed as in [13j are expected to account for the factoriz- 
able electroweak corrections, using the effective Born approximation [23], [21]. The non- 
factorizable contribution, very relevant in the large invariant mass region (Q = 160 GeV 
and above) are estimated to be much larger [23] . 

In all our studies we have fixed the energy of the collision to be = 14 TeV. 

The hadronic structure function Wy{t, Q"^) is represented by a convolution product be- 
tween the parton luminosities $]^(a;, and the Wilson coefficients Ajj(a;, Q^, /i^) 
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where the luminosities are given by 



X 



<^ij{x,fl%fil) = 

and the Wilson coefficients depend from both scales 

oo 



[fi ® fj] {x, /iL /^f) 



^2) 



n=0 



The explicit expressions of the hard scatterings coefficients have been taken from [12j and 
implemented in Candia. Moving to the parton densities, these are decomposed into their 
singlet (S) and nonsinglet (NS) contributions starting from the explicit expression 
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and after an expansion, one identifies, as usual, the convolution products S®S, NS®NS 
and S ® NS. 

As we have already mentioned, in each of this sectors we are entitled to implement 
evolved pdf 's of different accuracy, according to the classification presented in the previous 
section. Summarizing, we have, for the nonsinglet sector: 1) exact solutions of the NNLO 
exact equation; 2) exact solution of the NNLO truncated equation; 3) truncated solution 
of the NNLO truncated equation, while for the singlet case we have only the option of 
the K'-truncated solutions. 

As we have already explained, we work with the equations written according to form-1, 
which has a single expansion parameter (k'). This implies that the parton luminosities can 
be of a varying accuracy depending on the type of the solutions. The numerical analysis 
of these choices is very involved for realistic distributions, as we are going to discuss next. 
We remark that there are differences between the iterated solutions of type-1 and the 
hrute force solutions or the exact solutions, which are also available in the nonsinglet 
case. We have tried to answer this subtle point by showing in Fig. [1] the results for the 
cross sections determined at NLO and at NNLO using as input the MRST conditions 
taken from the grid, evolved by us using different sets of solutions. We recall that the 
initial condition /ig = 1 GeV^ means that we are using the MRST parametric input [IQJ. 
We have defined the "asymptotic solution" to be aasymp, built using the exact solution in 
the nonsinglet sector and a truncated solution in the singlet, with the index of truncation 
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K sufficiently large so that an asymptotic value for the logarithmic expansion (k = 8) is 
obtained. We plot the percentage difference, normalized as shown in the figure, between 
truncated solutions of a varying k index and this asymptotic cross section. It is clear, 
from this analysis, that the iterates of fixed accuracy, expanded in powers of \og{as)i do 
not converge to the asymptotic solution but give cross sections that differ by a small but 
finite amount from that. This is quantified to be of the order of 0.1 — 0.5% at the energy 
reported in the plot. This estimate is subject to change as we vary the energy scale and 
the model of the initial conditions. On the basis of this result, we may reasonably assume 
that the sequence of truncations, respect to the hrute force solution, or exact solution, 
should be of the order of a percent or so. This could be quantified better using a numerical 
code that solves the DGLAP by direct discretization, which is not available to us. From 
this point on, all the analysis that follows is going to be based on the implementation 
of cTasymp, as defined above. More details concerning the difference between truncated 
and asymptotic solutions, a critical analysis of these results and of their implications for 
precision studies of the parton model at NNLO will be presented below. 

5 Numerical Analysis: Comparison with the Les Houches 
and the MRST Models 

We start presenting in this section our comparisons between the results for the evolution 
and the cross sections obtained using Candia against those of other implementations. 
In doing this we have made sure that the same conditions are kept in regard both to the 
treatment of the heavy flavors and of the initial conditions when running the different 
codes. In particular, the parameters of the runs have been selected so to generate either 
truncated solutions or asymptotic solutions, as specified above. 

5.1 Comparisons with the Les Houches benchmarks 

We start our comparison using as initial conditions those presented in the Les Houches 
Model [21], which have been used to determine some benchmarks for the evolution. The 
implementations that we compare, in this case, are those of Candia and Pegasus, the 
latter running with the option IMODEV= 1. This option generates exact solutions of 
the evolution equations by using a large sequence of truncates in Mellin space, with both 
parameters k and k' large, according to the [/-ansatz (12^ . The heavy quarks have been 
treated according to the VFN scheme. 

In the Les Houches model pT] the input distributions mimic the CTEQ5M [25] pa- 



29 



0.0145 



0.014 



0.0125 



0.012 



NLO xu„(x), = Hp= 100 GeV, Nj= 4 
K=l,2,3,4,5,6,7-truncated xu„(x) 




0.0001 0.000105 0.00011 0.000115 

X 

(a) xuv{x, jjLp) at NLO 



NLO xg(x), = Hp= 100 GeV, Np 4 
K=l,2,3,4,5,6,7-truncated xg(x) 
LO 



82 
0.0001 



0.000105 



0.00011 



0.000115 



(c) xg{x,^F) at NLO 



0.0155 



0.015 



0.0145 - 



0.014 



0.0135 



0.013 



0.00012 




0.00012 



83 
0.0001 



NNLO xu^(x), = Hp= 100 GeV, Np 4 
K=2,3,4,5,6,7-truncated xu„(x) 



0.0001 0.000105 0.00011 0.000115 

X 

(b) xUy{x,iJ,F) at NNLO 



NNLO xg(x), = Hp= 100 GeV, Np 4 
K=2,3,4,5,6,7-truncated xg(x) 



0.000105 



0.00011 



0.000115 



(d) xg{x,HF) at NNLO 



Figure 2: Asymptotic and truncated pdf's for the valence up-quark and for the gluons 
at NLO and NNLO with nj^ — np — 100 GeV. We have selected the Les Houches input 
and a fixed number of flavors, Nf — 4. The small range for x has been chosen to resolve 
among the various predictions. 
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rameterization and are used regardless of the order of the evolution equations. They are 
given by 

xu,,{x) = 5.107200x°-^(l - x)^ 
xd^{x) = 3.064320a;°-^(l - x)^ 
xg{x) = 1.700000x-°'^(l - xf 
xd{x) = 0.1939875x"°-i(l - xf 
xu{x) = (1 — x)xd{x) 

xs{x) = xs{x) = 0.2x{u + d){x) , (84) 

and the running coupling has the value as(/^|;o ~ 2GeV^) = 0.35. Our implementation 
in Candia of the heavy thresholds, in this case, follows exactly the one described in 
[19]. To show the very good agreement between our method of solution and Pegasus 
we detail the results for all the sectors. We have included both the numerical values for 
the pdf's and the LO, NLO and NNLO predictions for the cross sections obtained by the 
two different implementations of the evolution. Tables [1]- [6] show the gluon and u-quark 
distributions using the two evolutions at the various orders. In both cases we keep the 
"asymptotic" mode (IM0DEV=1 for Pegasus) and the asymptotic solutions in Candia, 
with the nonsinglet treated using the exact iterated ansatz. It can be noticed that the 
differences are very small for all the densities up to NNLO. They can be read directly from 
the Tables ( fT|3f5i) since x6f{x) are the relative differences normalized to the Pegasus 

determination, i.e. x6f{x) = (a;/(x)pE(.^sus - ^/(^)Candia)/^/(^)pegasus- ^he 
percentage differences for the gluon densities are 0.2 % or smaller at NLO, 0.4% and 
smaller at NNLO. In the kinematical region relevant for the LHC they stay around 0.1% 
at NNLO. The valence u-quark distributions, at NNLO, reach at most 1% at a; = 10~^, 
while they are about 0.1 — 0.2% at a; = 10~^. Coming to the cross sections, the differences 
between the two determinations are pretty small. They essentially coincide at LO, they 
are about 0.6% at NLO, while they are about 0.3% and below at NNLO (see Tables 
andE]). 

5.2 Truncations and asymptotic solutions 

The reader can find in a sequence of 8 tables (see Tabs. I21II22|23|24|I25|26II27II28P added at 
the end of this work detailed numerical results for the various truncated solutions and for 
the corresponding asymptotic solution in the Les Houches model and in a realistic model, 
MRST [lU]. 

We show in Figs. ([2]) four plots of the valence up-quark and of the gluon distributions 
for various k values. The small range of variability in x has been chosen so to render the 
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differences in the plots visible, since they are quite small. The various k solutions converge 
toward the asymptotic solution as the index of truncation increases. We show the exact 
(for the valence up-quark distribution) or the asymptotic solution (for the gluon density) 
and the various truncated solutions for several k values. In the case of the Les Houches 
model, table EH and [22] are particularly significant, since these show for the nonsinglet 
the existence of a difference between the exact solutions, that performs a resummation 
of the log(as), and the sequence of truncated solutions, which reach saturation at k = 6. 
The differences for the valence up-quark distribution {xUv{x)) at NLO vary from 1% at 
{x = 10~^) to 0.7% {x = 10"'^), growing larger at x = 10~^, where they reach 7%. 
This last value is presented only for comparison, although it is not relevant at the LHC 
Moving to NNLO, the differences are about 4% at x = 10~^, 3% at x = 10""^, decreasing 
to 0.4% at X = 10~^. They become significant at large x values, being around 9% at 
X = 0.9. These determinations, of course, need to be tested in the related cross sections 
in order to appreciate their real impact. As we have already shown in Fig. [T]the various 
determinations stay below 1% for Q = Mz- Even if these differences are not big, they will 
become more significant as the determination of the pdf 's is going to improve in the near 
future, using the large amount of data coming from the LHC. This will allow to reduce 
the errors on the pdf 's and, therefore, on the cross sections. As we are going to show next, 
these errors remain, at the moment, larger than the theoretical indetermination coming 
from the choice of the solution, at least in the region that we have explored. In the gluon 
sector (see Tabl2Sl the situation seems to improve, and the differences stay below 1 
% in all the x-range, but this can be misleading: asymptotic and truncated solutions in 
the singlet sector are in fact both determined by the same logarithmic ansatz. 

A similar analysis has been performed for the MRST model. In this case we perform 
the evolution using Candia, the MRST input and a treatment of the heavy flavors exactly 
as in MRST, with the thresholds for the heavy quarks chosen as in [10]. Also in this case 
truncated solutions and asymptotic solutions show a small difference, both for the valence 
distributions and for the singlet ones. We show in tables [2^26|27l and [2H1 results for the 
various K-truncated (accurate) solutions. 

For instance, in the case of the gluon density, if we choose k = 3 (3rd truncated 
solution), at X = 10"'^ the difference in the gluon density respect to the asymptotic 
solution is about 0.01% at NNLO, which appears to be small, but can easily grow to 0.5% 
or so if would let a brute force solution replace the asymptotic determination. In fact the 
valence u-quark distribution, whose asymptotic value is supposed to be pretty close to 
the exact value, shows more substantial differences. For instance, at NLO, for x = 10~^ 
the same truncated ansatz {k = 3) differs from the exact one by 2.6%. At NNLO in the 
more relevant region of x (0.01-0.1) is about 2% and below. The differences grow bigger 
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Candia vs Pegasus PDFs at LO, Les Houches input, VFN scheme, Q ^ ^ i^R = 100 GeV 
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Table 1: Comparison between the pdf's obtained using Candia versus those ob- 
tained using Pegasus and the normahzed differences, ex.: Sxg{x) = — 
x(7(x)PegasuS|/^^(^)Pegasus LO. 



at larger x-values, for instance they are 9% for x = 0.5 at NNLO. 

Coming to the cross sections obtained by the various truncated solutions, these are 
shown in two tables (see Tabs. [ ^[HUj) . which summarize these studies at NLO and NNLO 
respectively. Using again the k = 3 solution, for Q = Mz the NLO determination differs 
by 0.2% compared to the asymptotic one. They tend to grow at larger Q-values, 0.4% at 
Q = 200 GeV (NNLO). 

There are some conclusions that we can draw from this analysis. We clearly have 
several ways to choose the solution and by doing so we make errors which are around 
1%. They tend to grow as Q increases, at larger invariant mass of the lepton pair, where 
we get more sensitive to larger x-values. This theoretical errors may grow slightly bigger 
at very large Q-values, say for Q around 1 TeV or so, where we need specific studies of 
that kinematical region, since we could expect that extra neutral interactions be found. 
It is important, however, to remind that the DY cross section is anyhow quite sensitive 
to the behavior of the hard scatterings around x = 1, as pointed out in [12]. This implies 
that various determinations may differ already at percent level because of the different 
treatment of the edge-point region in the Bjorken variable even for moderate Q values. 
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pb/GeV]. Candia vs Pegasus with Les Houches input. 
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Table 2: Comparison between the cross sections obtained using Candia and Pegasus 
at LO. 



Candia vs Pegasus PDFs at NLO, Les Houches input, VFN scheme, Q = [if = \i-R = 100 GeV 
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Table 3: Pdf 's obtained in the two evolutions at NLO 
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do-AT/^o/dQ [pb/GeV]. Candia vs Pegasus with Les Houches input. 
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Table 4: NLO cross sections obtained using Candia and Pegasus using the Les Houches 
model. 



Candia vs Pegasus PDFs at NNLO, Les Houches input, VFN scheme, Q = hf = IJ-r = 100 GeV 


X 


„„/„\(jANDIA 


a;5(a;)^EGASUS 


6xg{x) 


Candia 


a;«,(a;)^EGASUS 


6xUy {x) 


le-05 


2.1922-10+2 


2.2012 - 10+2 


4.1108-10-3 


3.0823 - 10-3 


3.1907-10-3 


3.3962-10-2 


0.0001 


8.8486 - 10+1 


8.8804- 10+1 


3.5856-10-3 


1.3871 - 10-2 


1.4023- 10-2 


1.0811 - 10-2 


0.001 


3.0319-10+1 


3.0404 • 10+1 


2.8106-10-3 


6.0060 - 10-2 


6.0019-10-2 


6.9117- 10-4 


0.01 


7.7785 - 10+0 


7.7912- 10+0 


1.6326-10-3 


2.3287-10-1 


2.3244- 10-1 


1.8584- 10-3 


0.1 


8.5284 - 10-1 


8.5266 - 10-1 


2.1595 - 10-4 


5.4977 - 10-1 


5.4993 - 10-1 


2.9526 - 10-4 


0.2 


2.4183-10-1 


2.4161 - 10-1 


9.1195-10-4 


4.8313 - 10-1 


4.8323 - 10-1 


2.0148-10-4 


0.3 


7.9005 - 10-2 


7.8898 • 10-2 


1.3515 - 10-3 


3.4629-10-1 


3.4622- 10-1 


1.9857-10-4 


0.4 


2.5636 - 10-2 


2.5594-10-2 


1.6452 • 10-3 


2.1711 - IQ-i 


2.1696- 10-1 


6.7488 - 10-4 


0.5 


7.6538 - 10-3 


7.6398 - 10-3 


1.8314-10-3 


1.1883-10-1 


1.1868-10-1 


1.2434 - 10-3 


0.6 


1.9439-10-3 


1.9401 - 10-3 


1.9844 - 10-3 


5.4753 - 10-2 


5.4652 - 10-2 


1.8520-10-3 


0.7 


3.7162 - 10-4 


3.7080 - 10-4 


2.2059 - 10-3 


1.9537-10-2 


1.9486 - 10-2 


2.6105-10-3 


0.8 


4.1248-10-^ 


4.1141 - 10-^ 


2.5990 - 10-3 


4.4306 - 10-3 


4.4148 - 10-3 


3.5750 - 10-3 


0.9 


1.1766- lO-'^ 


1.1722-10-6 


3.7723 - 10-3 


3.3696 - 10-4 


3.3522 - 10-4 


5.1816-10-3 



Table 5: NNLO pdf's determined with Candia and Pegasus using the Les Houches 
model. 
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dcTATATLo/dQ [pb/GeV]. Candia vs Pegasus with Les Houches input. 
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3.6106 - 10-3 


70.0938 


4.4033 ■ 10+° 


4.4184 • 10+° 


3.4110 - 10-3 


80.1407 


8.9241 ■ 10+° 


8.9527-10+° 


3.1936-10-3 


90.1876 


3.3570 ■ 10+2 


3.3669 ■ 10+2 


2.9388 - 10-3 


91.1876 


5.4905 ■ 10+2 


5.5067-10+2 


2.9299 - 10-3 


92.1876 


3.3344 ■ 10+2 


3.3441 - 10+2 


2.8919 - 10-3 


120.0701 


1.0249 ■ 10+° 


1.0274-10+° 


2.4285 - 10-3 


146.0938 


2.8527-10-1 


2.8590-10-1 


2.1826 - 10-3 


172.1175 


1.2295 ■ 10-1 


1.2319- 10-1 


1.9887- 10-3 


200.0000 


6.0923 ■ 10-2 


6.1029- 10-2 


1.7369- 10-3 



Table 6: NNLO cross sections in the two evolution methods. 



6 Other comparisons with the MRST evolution 

Now we perform a comparison between the various cross sections obtained using the pdf 's 
evolved by MRST and the same distributions, taken at their starting value, but evolved by 
us using Candia. These studies are performed using in Candia the asymptotic solutions 
in the singlet and non singlet sectors. We use the MRST input in a grid form with an 
initial scale /Xq = 1.25 GeV2, ^/S = 14 TeV and with = /i^ = Q'^. The choices for 
the thresholds of the heavy flavors have been chosen in Candia to coincide with those 
reported by MRST. For this reason we have used for comparison the asymptotic solution 
and the VFN scheme. The relative variations are computed respect to the MRST value 
and are indicated in the columns labeled as 6(J. Also in this case we show in 3 tables 
([7}in]) the results for the LO, NLO and NNLO cross sections. The differences between our 
prediction and the MRST result for the total cross sections are around 1 per cent or below 
at LO, vary from 0.02% to 0.3% at NLO and are 2.6% and below at NNLO. In this case 
the maximum difference has been found for Q = 50 GeV. These differences, clearly, affect 
the values of the ii'-factors, as we are going to discuss below, which in our evolution are 
larger compared to those of MRST. 

We perform some more tests using Vrap [13] and compare the results against those 
of CandiAdy for the calculations of the hard scattering piece. In the results given below 
^Vrap I'sfers to the hard scatterings for the invariant mass distributions computed using 
Vrap, while ^MRSTrt refers to the luminosities using one of the MRST inputs 
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evolved using Candia. Similarly, ^mrst^^oi denotes the luminosities predicted by MRST 
with their evolution. In this case the original scale is not indicated, but the grid scale 
Ho — 1.25 GeV^ is the first available point, at which the evolution with Candia is also 
interfaced. 

The NLO total cross sections in [pb/GeV] at the peak Q = Mz and at = 14 TeV 
using the MRST inputs with jJ^l = I GeV^ and //q = 1.25 GeV^, in the various cases, are 
given by 



^y™>^Mii5T_, = 501.96, 



Krap ® ^MRST^^^^^^ = 505.87 from /.^ = 1 GeV^ , 
^vrS ® "^^^^^Candia = ^'^"^ = ^-^^ ^""^^ 



^ciNDlAw®*Mi.5T... = 501.72, 



^CANDlAzpy ® ^^^^^CANDIA = ^^^'^^ from /.^ = 1 GeV^ , 
^ciNDlAw ® ^^^^^Candia = ^°2.42 from /.g = 1.25 GeV^ (85) 
with differences that stay well below 1%, while at NNLO we obtain 



^y.T ® "^MRST^.^, = 490.51 , 

^vZi"" ® *^«5^Candia ^ ^^^-^^ ^'""^ = ^ ' 
'^yi^'' ® ^^^^^Candia ^ ^^^-^^ ^™ = 1-25 GeV^ , 
^CrNDiAw®^Mi.5T...= 488.22, 

-Candiaw ^ ^^^^^Candia = ^^^-^l ^^^^ ^0 = 1 ' 

^cInDIAW ® ^^^^^CANDIA = ^^""^ ^0 = 1-25 ^S^) 

cross sections that differ approximately by 2%. The reduction of the cross section in 
Candia is more remarked compared to MRST and is due to the evolution. 



6.1 The renormalizat ion/factorization scale dependence of the 
cross section 

An interesting aspect of the prediction of the QCD observables is their factorization and 
renormalization scale dependence. We will denote hy /ip and /ir the two scales. The 
dependence is important and appears both in the hard scatterings and in the evolved 
pdf's, using the modified NNLO kernels defined above. The optimal choices for these 
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do-Lo/dQ 


pb/GeV]. Candia vs MRST evol. with MRST input, /xg = 1.25 GeV^ 


Q [GeV] 


^(Jandia 


^MRST 




50.0000 


5.6629 ■ 10+° 


5.7110 ■ 10+° 


8.4230 ■ 10-3 


60.0469 


3.4301 ■ 10+° 


3.4692 ■ 10+° 


1.1274-10-2 


70.0938 


3.1248 ■ 10+° 


3.1646 • 10+° 


1.2583 ■ 10-2 


80.1407 


6.4675 • 10+° 


6.5540 • 10+° 


1.3191 ■ 10-2 


90.1876 


2.4859 • 10+2 


2.5189 • 10+2 


1.3086 • 10-2 


91.1876 


4.0723 • 10+2 


4.1261 • 10+2 


1.3059 • 10-2 


92.1876 


2.4767 • 10+2 


2.5094 • 10+2 


1.3033 • 10-2 


120.0701 


7.6837 • 10-1 


7.7755 • 10-1 


1.1796 • 10-2 


146.0938 


2.1196-10-1 


2.1415 • 10-1 


1.0240 • 10-2 


172.1175 


9.0345 • 10-2 


9.1149 • 10-2 


8.8207 • 10-3 


200.0000 


4.4185 • 10-2 


4.4504 • 10-2 


7.1679 • 10-3 



Table 7: LO cross section for Drell-Yan obtained by Candia using the MRST input and 
the evolved MRST pdf 's 



dcTATio/dQ [pb/GeV]. Candia vs MRST evol. with MRST input, nl = 1.25 GeV2 


Q [GeV] 


^Candia 

'^NLO 


^MRST 
'^NLO 


Sctnlo 


50.0000 


6.8119 • 10+° 


6.8100 • 10+° 


2.7680 • 10-^ 


60.0469 


4.1552 • 10+° 


4.1521 • 10+° 


7.5793 • 10-^ 


70.0938 


3.8110 • 10+° 


3.8080 • 10+° 


8.1120-10-^ 


80.1407 


7.9371 • 10+° 


7.9287 • 10+° 


1.0526 • 10-3 


90.1876 


3.0657 • 10+2 


3.0615 • 10+2 


1.3656 • 10-3 


91.1876 


5.0242 • 10+2 


5.0172 • 10+2 


1.3903 • 10-3 


92.1876 


3.0569 • 10+2 


3.0526 • 10+2 


1.4133 • 10-3 


120.0701 


9.5677 • 10-1 


9.5496 • 10-1 


1.8964 • 10-3 


146.0938 


2.6562 • 10-1 


2.6504 • 10-1 


2.1997 • 10-3 


172.1175 


1.1382 ■ 10-1 


1.1356-10-1 


2.2278 • 10-3 


200.0000 


5.5940 • 10-2 


5.5778 • 10-2 


2.9044 • 10-3 



Table 8: NLO cross section for Drell-Yan obtained by Candia using the MRST input 
and the evolved MRST pdf 's 
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Figure 3: 
MRST. 



Various X-factors obtained with the evolution performed by Candia and 
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^(yNNLo/^Q [pb/GeV]. Candia vs MRST evolution with MRST input, /ig = 1.25 GeV^ 


Q [GeV] 


^(Jandia 

^NNLO 


-.MRST 
^NNLO 




50.0000 


6.4935 • 10+° 


6.6707-10+° 


2.6560- 10-2 


60.0469 


3.9997-10+° 


4.0961 - 10+° 


2.3534-10-2 


70.0938 


3.6962 • 10+° 


3.7743 - 10+° 


2.0678-10-2 


80.1407 


7.6755 ■ 10+° 


7.8198 - 10+° 


1.8455-10-2 


90.1876 


2.9325- 10+2 


2.9827-10+2 


1.6834-10-2 


91.1876 


4.8006 ■ 10+2 


4.8822 - 10+2 


1.6702-10-2 


92.1876 


2.9179-10+2 


2.9671 - 10+2 


1.6575-10-2 


120.0701 


9.0411 - 10-1 


9.1687-10-1 


1.3918-10-2 


146.0938 


2.5267- 10-1 


2.5567-10-1 


1.1714-10-2 


172.1175 


1.0938-10-1 


1.1049 - 10-1 


1.0028-10-2 


200.0000 


5.4431 - 10-2 


5.4876 - 10-2 


8.1092-10-3 



Table 9: NNLO cross section for Drell-Yan obtained by Candia using the MRST input 
and the evolved MRST pdf's 



scales are identified in the region of stability of the cross section, usually a small plateau in 
a multi-parameter space, which can be searched numerically. In the case of the Higgs, for 
instance, a rather general analysis of the structure of these surfaces for specific observables 
(such as the total cross section for Higgs production and the corresponding i^"-factors) 
through NNLO has been given. There one can show, but the result is quite general, that 
the concavity of the bidimensional surfaces describing the cross sections, plotted in terms 
of the two independent scales, changes sign as we move from leading to next-to-leading 
order [26]. We show in Figs. H] and [5] global plots of the DY cross section near the Z peak 
for the two models Alekhin and MRST, evolved by the same authors, and zooms of the 
peak region (FigJ^l), in which the LO, NLO and NNLO contributions are resolved in great 
detail. Here we have set the factorization scale to be Q (/i^ = Q). In two following plots. 
Fig. [6] and [TJ we show instead the variation of the same cross section using an evolution 
provided by Candia at LO, NLO and NNLO of the MRST input from the grids {^l = 1.25 
GeV2)and we have varied /x^ and fi^. As we have already mentioned, our analysis includes 
all the fiR dependence (see Tabs. [TU|llfl2lll3lll4p . coming both from the pdf's and from the 
hard scatterings. The first is usually not reported in the standard parameterizations such 
as MRST and Alekhin. The variation has been performed setting, for each (fixed) value 
of Q, fip = Q and studying the variation of the renormalization scale in the ratio kp = 
fj^ji/ fj^'p, which has been taken to vary between 1/2 and 2. The decreased dependence of 
the result on the spurious scales of the process as we move toward the NNLO predictions 
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from the LO ones are quite visible. This is particularly easy to see from Fig. [71 The 
options fiR = hf and Q ^ Hf are shown both for Alekhin's and the MRST inputs (Figs. [8] 
and where a zoom of the region above the Z peak and of the tail of the cross section 
are presented. The region covered is quite small (100-110 GeV) so to allow to discern 
between the various results. The bands of variations of the LO, NLO and NNLO results 
can be identified by a close look at these figures. One can see immediately the reduced 
sizes of these bands as we increase the perturbative order of accuracy. The same bands 
are shown right on the peak of the Z in Fig. One can immediately notice that the 
the NNLO variations take place right inside the NLO error band for the Alekhin model 
(Fig. [9] (a)), while they overlap at the edge in the MRST model (Fig. [9] (b)). Regarding 
the precise size of these variations, these can be inferred from the corresponding tables. 
The range explored in our analysis (1/2 < fc^ < 2) is somehow smaller than that explored 
in [13], but includes the entire dependence on the renormalization scale of the pdf's. Being 
the evolution rather important in the determination of the NNLO total cross section, it 
is clear that also the ^ir dependence on the evolution is not negligible. The two cases 
fJ'R < fJ'F and fiji> fiF are characterized by substantially different excursions in range. In 
the first case the variations, at LO, are from 25 % at 50 GeV down to 10 % for Q = 200 
GeV, while for kF = 2 they are more moderate (from 17 % down to 8 %). At NLO the 
excursions are approximately from 11% down to 6% in the same range of Q, for both cases 
of kF- The variations at NNLO can be found in[T21 and are in the range of 1-3 %. We 
have also shown in tab. [T3] and [T3] results for the scale dependence when we remove 
in the pdf's, by equating to fiF and keep them separate only in the hard scatterings. 
The range of variation are sensibly reduced especially at lower values of Q, with a drastic 
reduction especially around the peak. The reduction in the variation is by a factor of 
10 less: from about 10% down to less than 1%. On the peak the NNLO variations are 
between 0.1 and 0.03 %. It is clear from this results that the fiR scale dependence coming 
from the evolution is pretty relevant and, in a complete analysis of the stability of the 
NNLO corrections can't be forgotten. 

6.2 The K factors 

We have summarized in Fig. [3] four plots of the behavior of the 3 i^-factors K = 
o'NNLo/o'NLO, Ki = aNLo/o'LO and K2 = aNNLo/c^LO obtained using Candia and the 
MRST evolution. These are shown as a function of Q, and evaluated at the center of mass 
energy of = 14 TeV. The dependence of the results on the evolution is significant. In 
fact, from Fig. [3] it is evident that while the shapes of the plots of the i^-factors are simi- 
lar, there are variations of the order 2%, in the results using the two different evolutions. 
Both in the evolution performed with Candia and in the MRST evolution we use the 
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daLo{Q,fJ'F,fJ'R)/dQ [pb/GeV]. Candia evolution with MRST input, fj,^ = 1.25 GeV^ 


Q [GeV] 




kp = 1 


kp = 1/2 






50.0000 


6.0789 • 10+° 


5.6629 • 10+° 


7.0877 • 10+° 


2.5159- 10+1 


1.7291 • 10+1 


60.0469 


2.8718 • 10+° 


3.4301 - 10+° 


4.2241 - 10+° 


2.3148- 10+1 


1.6276-10+1 


70.0938 


2.6438 • 10+° 


3.1248-10+° 


3.7952 - 10+° 


2.1454- 10+1 


1.5394-10+1 


80.1407 


5.5230-10+° 


6.4675 - 10+° 


7.7598 - 10+° 


1.9981- 10+1 


1.4604-10+1 


90.1876 


2.1411 • 10+2 


2.4859 • 10+2 


2.9498 • 10+2 


1.8662 • 10+1 


1.3869 • 10+1 


91.1876 


3.5104 • 10+2 


4.0723 • 10+2 


4.8272 • 10+2 


1.8538-10+1 


1.3798 • 10+1 


120.0701 


6.7639-10-1 


7.6837-10-1 


8.8706-10-1 


1.5446 - 10+1 


1.1972-10+1 


146.0938 


1.8945-10-1 


2.1196 - 10-1 


2.4008 - 10-1 


1.3267- 10+1 


1.0620-10+1 


172.1175 


8.1799-10-2 


9.0345-10-2 


1.0070-10-1 


1.1463- 10+1 


9.4593 - 10+° 


200.0000 


4.0486 • 10-2 


4.4185-10-2 


4.8524 • 10-2 


9.8201 • 10+° 


8.3716 • 10+° 



Table 10: Study of the variation of the LO cross sections with respect to kp = n\/ 
Here we choose Q — /if and the 11%/ variation is also included in the pdf's evolved 
with Candia. 



dcrNLo{Q,fJ'F,m)/dQ [pb/GeV]. Candia cvol. with MRST input, /x§ = 1.25 GeV2 


Q [GcV] 


'T^^^CQ) kp = 2 


fT^^^CQ) kp = 1 


a^'-'^iQ) kp - 1/2 






50.0000 


6.0789 - 10+° 


6.8121 - 10+° 


7.5694 - 10+° 


1.1116-10+1 


1.0763 - 10+1 


60.0469 


3.7343 - 10+° 


4.1554-10+° 


4.5906 - 10+° 


1.0473-10+1 


1.0134-10+1 


70.0938 


3.4443 • 10+° 


3.8112 • 10+° 


4.1920 • 10+° 


9.9929 • 10+° 


9.6262 • 10+° 


80.1407 


7.2077 • 10+° 


7.9374 • 10+° 


8.6992 • 10+° 


9.5986 • 10+° 


9.1922 • 10+° 


90.1876 


2.7965 - 10+2 


3.0658 - 10+2 


3.3480 - 10+2 


9.2047 - 10+° 


8.7856 - 10+° 


91.1876 


4.5849 - 10+2 


5.0243 - 10+2 


5.4848 - 10+2 


9.1649-10+° 


8.7461 - 10+° 


120.0701 


8.8285-10-1 


9.5681-10-1 


1.0345-10+° 


8.1197-10+° 


7.7303 - 10+° 


146.0938 


2.4702 - 10-1 


2.6563-10-1 


2.8526-10-1 


7.3907-10+° 


7.0048 • 10+° 


172.1175 


1.0G54 ■ 10-1 


l.i:-!<S2 ■ i()-i 


1.21.5.5 ■ 10-1 


().7!)40 ■ 1U+" 


().;]<)78 ■ 10-" 


200.0000 


5.2673 • 10-2 


5.5942 • 10-2 


5.9441 • 10-2 


6.2547 • 10+° 


5.8436 • 10+° 



Table 11: Study of the variation of the NLO cross sections with respect to kp = iJi\/ iJ?p. 
Here we choose Q = i^f and the fJ^ji/ fJ^p variation is also included in the pdf's evolved 
with Candia. 
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<iaNNLo{Q,IJ'F,IJ'R)/dQ [pb/GeV]. Candia evolution with MRST input, i^l = 1.25 GeV^ 


Q [GeV] 






^NNLOf^Q^ Up = 1/2 






50.0000 


6.2855 • 10+" 


6.4940 • 10+" 


6.5465 - 10+" 


8.0790-10-1 


3.2107- 10+" 


60.0469 


3.8626 • 10+" 


3.9989 - 10+" 


4.0503 - 10+" 


1.2850- 10+" 


3.4090 - 10+" 


70.0938 


3.5635 • 10+" 


3.6948 - 10+" 


3.7557 - 10+" 


1.6481 - 10+" 


3.5528 - 10+" 


80.1407 


7.4312 • 10+" 


7.6740 • 10+" 


7.7729 - 10+" 


1.2890- 10+" 


3.1640- 10+" 


90.1876 


2.8672 • 10+2 


2.9335 • 10+2 


2.9409 • 10+2 


2.5092-10-1 


2.2615 • 10+" 


91.1876 


4.6986 • 10+2 


4.8027 • 10+2 


4.8094 • 10+2 


1.4037-10-1 


2.1663 • 10+" 


120.0701 


9.0247- 10-1 


9.0552-10-1 


8.8841 - 10-1 


-1.8901 - 10+" 


3.3748 - 10-1 


146.0938 


2.5416- 10-1 


2.5318-10-1 


2.4664- 10-1 


-2.5808- 10+" 


-3.8826- 10-1 


172.1175 


1.1060- 10-1 


1.0963-10-1 


1.0633-10-1 


-3.0129- 10+" 


-8.8480-10-1 


200.0000 


5.5290-10-2 


5.4572-10-2 


5.2727- 10-2 


-3.3809- 10+" 


-1.3157- 10+" 



Table 12: Study of the variation of the NNLO cross sections with respect to kp = 

Here we choose Q — fxp and the iiji/ /J^p variation is also included in the pdf's evolved 

with Candia. 



dcTNLoiQ, fJ-F, fJ'R)/dQ [pb/GcV]. Candia evolution with MRST input, /Xq = 1.25 GcV2 aikp) <S> ^(fJ-p) 


Q [GcV] 


fT^^^CQ) kp = 2 


cj^^^iQ) kp = 1 


a^'-^iQ) kp = 1/2 




SakF=2% 


50.0000 


6.7636 - 10+" 


6.8121 - 10+° 


6.8667- 10+" 


8.0201-10-1 


7.1156-10-1 


60.0469 


4.1271 - 10+" 


4.1554- 10+" 


4.1871 - 10+" 


7.6402 - 10-1 


6.8044-10-1 


70.0938 


3.78()3 ■ 10+" 


3.8112 ■ 10+" 


3.8390 ■ 10+" 


7.3124 ■ IQ-i 


().r>:m ■ 10-1 


80.1407 


7.8873 - 10+° 


7.9374 - 10+° 


7.9933 - 10+" 


7.0434-10-1 


6.3080-10-1 


90.1876 


3.0469 - 10+2 


3.0658 - 10+2 


3.0869 - 10+2 


6.8790-10-1 


6.1737-10-1 


91.1876 


4.9934 - 10+2 


5.0243 - 10+2 


5.0589-10+2 


6.8677-10-1 


6.1649-10-1 


120.0701 


9.5104 - 10-1 


9.5681 - 10-1 


9.6321 - 10-1 


6.6868- 10-1 


6.0294-10-1 


146.0938 


2.6405 • 10-1 


2.6563-10-1 


2.6738 • 10-1 


6.5806 • 10-1 


5.9519-10-1 


172.1175 


1.1315-10-1 


1.1382-10-1 


1.1456-10-1 


6.5014-10-1 


5.8952 • 10-1 


200.0000 


5.5615 • 10-2 


5.5942 • 10-2 


5.6303 • 10-2 


6.4531 • 10-1 


5.8453 • 10-1 



Table 13: Study of the variation of the NLO cross sections with respect to kp — 

Here we choose Q — /ip and the IJ>%/ /ijr variation is only included in the hard scattering 

piece. 
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dcrNNLo{Q,IJ'F,IJ'R)/<iQ [pb/GeV]. Candia evolution with MRST input, /Xg = 1.25 GeV^ ^{kp) ^{i^f) 










da, 1 




K(\ nnnn 
ou.uuuu 




a AQAC\ . 1 n+0 


o.oooy • lu 


D.yioo • lu 


R OKA A . 1 0~1 


fin DzLfiQ 




^ QQRQ . 1 0+^ 


zi oziRzi . 1 


1 9*^71 . 1 0+0 


1 nfizis . 1 0+0 


( u.uyoo 


Q fi/l /1 9 . 1 


Q AQAS 1 n+'^ 
o.uy^o • lu 


Q 7KAA 1 


1 fil 9"^ . 1 O+O 


1 '?7nA . 1 0+0 
1.0/ - lU 




7 ^Qfi'^ . 1 n+0 


7 f^7zlf^ . 1 n+0 




1 1 71 1 . 1 O+O 


1 01 99 . 1 0+0 




9 0*^1 t=l . 1 n+2 
z.yoio • ±u 






1 f\'^AA . 1 0~2 
X.DO^^ • lU 


a 70C7 . 1 n— 2 


Qi ia7« 


A 55nA9 . 1 n+2 


A 55097 . 1 n+2 


A 7Q77 . 1 0+2 


1 0'?SQ . 1 0~1 

— i.uooy • lu 


Q 074 1^ . 1 0~2 




Q 91 fifi . 1 n~i 


Q n'=»'^9 . 1 


R S^'^Q . 1 0~1 
o.oooy • lu 


9 99'^'^ . 1 0+0 
z.zziOO xu 


1 7S90 . 1 0+0 




9 ^QOQ . 1 n~i 


9 ^'^1 . 1 n~i 


9 zl'^SS . 1 

Zi. 4:000 XU 


— 9 SRI zi . 1 0+0 

^ . OO X4: xu 


-—9 '^'^fi7 . 1 0+0 


172.1175 


1.1256- 10-1 


1.0963-10-1 


1.0604- 10-1 


-3.2701 - 10+° 


-2.6699- 10+° 


200.0000 


5.6183-10-2 


5.4572 - 10-2 


5.2608 • 10-2 


-3.5989 - 10+° 


-2.9521 • 10+° 



Table 14: Study of the variation of the NNLO cross sections with respect to kp = fip- 
Here we choose Q — /if and the i-iji/ variation is only included in the hard scattering 
piece. 



same MRST input, choosing the initial scale //q = 1.25 GeV^, and the same treatment of 
the heavy flavors. On the Z resonance we get 

K{Mz) = {aNNLO ® ^ZlsT)/i^NLO ® ^Zr%) = 0-97 
K{Mz) = (a^NLO ® *crNDiA)/('^™ ® ^cTndia) = ^'^^ 

K{Mz) = {a^NLO ® ^AUkDK^NLO ® ^ILkhrn) = 0.98 (87) 

which corresponds to a reduction by 2.7% of the NNLO cross section compared to the 
NLO result, (MRST evolution) and larger for the Candia evolution, 4.4%, while for 
Alekhin is 1.5%. From the analysis of the errors on the pdf's to NNLO, for instance for 
the Alckhin's set, the differences among these determinations are still compatible, being 
the variations on the i^T-factors of the order of 4%. We will get back to this point in the 
next sections. Similar i^-factors can be introduced to study the variations from LO to 
NNLO. We obtain 

K^{Mz) = {aNNLO ® ^m1st)/{^lo ^mrst) = 1-18 
Ki{Mz) = {aNLO ® ^m''rst)/{^lo ® $mV) = 1-21 

MMz) = {a^Lo ® <^'Candia)/(^^o ® *g^,NDiA) = 1-23 , 

K^{Mz) = {CJNLO ® ^ILkhJIi^LO ® ^MekMn) = 1-23 , 
K2{Mz) = (a^iVLO ® *^;tLn)/(^LO ® ^TlekUn) = 1-21- 

corresponding to a growth around 17 — 23%. 
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6.3 The rapidity distributions 



Another cross section of relevance is the calculation of the rapidity distributions of the 
lepton pair in the final state at the resonance of the Z at NNLO. Since the number of events 
expected from Drell-Yan cross at LHC is large, the study of these distributions will be 
very important for partonometry. As we have already mentioned, the analysis presented 
here is going to be rather short, and we hope to return to this point in a separate work. 
We perform a numerical calculation of the differential cross section interfacing Vrap \13\ . 
which computes the hard scatterings, with Candia. At this point we recall that the 
rapidity of the vector boson Z is defined as 

y-hoJ^), (88) 



2 ^\E-p, 

where E and pz are respectively the energy and the longitudinal momentum of Z in the 
center of mass frame of the colliding hadrons. Integrating over this variable one obtains 
the total cross section as 

/•(l/2)lnl/r , 
J(l/2)lnr "-^ 



where 

dY 



f f 

^ = 5^/ / dXidX2fa'{Xi,Qyfll,fl%/fil)f^%X2,Qyfll,fil/fil) X 

xi,X2,Q'^/fiF,fin/fiF). (90) 



ab 

da 



dY 

Notice that the evolution implemented in Candia allows to analyze the renormaliza- 
tion/f actor ization scale dependence also in the evolution, which is not present in the 
MRST parameterizations. 

If we set the scales to be equal, np = hr and vary /i^? in the interval 1/2Q < /i_F < 2(5 
we obtain the results in Fig. fllOp . which differ from those obtained in [T3J by 2% due 
to the different implementation of the evolution. Using Candia and as initial condition 
the MRST grid input with /ig = 1.25 GeV^ the NNLO band and the NLO one are 
resolved separately. From Fig. ffTT]) it is clear that including the fJ^ji/ fJ^j? effects in the 
pdf's evolution, the dependence on fiji is quite sizeable at NLO, but is reduced at NNLO. 

We show in Fig. f|T2|) the plots of the variations of the rapidity distributions at the 
three orders and the corresponding pdf's errors for Alekhin's model and for MRST for 
Q = Mz- In both cases the reduction of the variation of the cross sections as we move 
toward higher orders is quite evident. We report also the errors on these distributions 
obtained in both models, which also get systematically smaller as the accuracy of the 
calculation increases. 
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200 



Q [GeV] 
(b) MRST Evolution 



Figure 4: Cross Sections in the region of the peak of the Z boson at LO, NLO, and NNLO 
obtained using the luminosities evolved respectively by Alekhin and MRST 
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Figure 5: Cross Sections in the region of the Z with a zoom in the peak region. 
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Q [GeV] 

(a) Candia LO evolution for MRST parametric input. 




Q [GeV] 

(b) Candia NLO evolution for MRST parametric input. 



Figure 6: Factorization vs Renormalization scale dependence of the cross section at LO, 
NLO with \/S = 14 TeV. The pdf 's have been evolved by using the MRST parametric 
input at nl — 1 GeV^ 
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50 60 70 80 90 100 110 120 

Q [GeV] 

Figure 7: Factorization vs renormalization scale dependence of the cross section at NNLO 
at \^ = 14 TeV. The pdf 's have been evolved by Candia using the MRST parametric 
input at /^o = 1 GeV^. 

7 The Cross Sections and the Errors 

We are now going to quantify the errors coming from the pdf 's on the differential cross 
section in the peak region of the Z by setting the condition /i/j = fip. The numerical 
determination of the errors is computationally very intensive. We perform the analysis at 
LO, NLO and NNLO for the case of the Alekhin's pdf's, and only at NLO for the case of 
the MRST model, since the error analysis in the latter case is not available at LO and at 
NNLO. We present our results in Figs. (|T3|) . f|T^ . f|T5l) and f|T6|) at typical LHC energy 
(v^= 14 TeV). 

Fig. [I3] (a) shows that the 2 error bands at NLO and NNLO intersect, though the 
average NNLO cross section is located outside the area covered by the error band in the 
fast fall off region. This feature of the result is shown more clearly in Figs. [TH and [TJl A 
zoom of the same region is shown in Fig. [161 

The calculation of the error bands has been done following the usual theory of the 
linear propagation of the errors. Starting from the errors on the pdf's known in the 
literature (see [ID] , [27] , [9] ) , we have generated different sets of cross sections. Then, the 
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(a) Zoom of the tail of the DY cross section above the peak of the Z using 
the evolution provided by Alekhin. 
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(b) As above, but for the MRST parameterizations, with input energy 
fil = L25 GeV2. 

Figure 8: Factorization scale dependence of the cross section at LO NLO and NNLO with 
a/S* = 14 TeV for two models. 
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(a) Alekhin's evolution. 
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(b) MRST evolution with input energy = 1.25 GeV^ 
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Figure 9: Factorization scale dependence of the cross section at LO NLO and NNLO with 
= 14 TeV. Zoom in the peak region. 
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Figure 10: Rapidity distributions obtained by changing 1/2Q < yU^ < 2Q. Here we choose 
Q = Mz = fip- The evolution is based on Candia using MRST grid input with /Iq = 1.25 
GeV^, while we used Vrap for the calculation of the hard scattering pieces 

error on the cross section has been calculated using the formula 



where ak is the k-th cross section belonging to a certain set, and N is the number of free 
parameters, which is 15 for MRST and 17 for Alekhin. 

We show in Table [15] the values for the cross section with values obtained by the 
best fits and the errors at the corresponding orders. It is evident that the relevance 
of the NNLO corrections is reduced at lower Q, given the actual quantification of the 
pdf's errors, since the NNLO corrections are not outside the error bands. The situation, 
however, changes beyond the resonance (120 GeV and above) , where it is clear that the 
cross section of best the fit at NNLO lays outside the error band, on the tail of the region 
that we analyze. The errors induced on the i^'-factors ( {K{Q) — 1)%), as one can easily 
figure out, are of the order of 4% on the peak {K{Mz) for the Alekhin set) from their 
best fit value, widening quite sharply that determination {K{Mz) = 0.98 ± 0.04). For 
Q = 50 GeV the percentile variation of the cross section in moving from NLO to NNLO is 
about 1.5 ± 4%. As a last example, for Q = 146 GeV we obtain a rate of variation of 3% 
(i^(146) = 0.97 ±0.03). Regarding the size of the errors at NNLO at various Q values, in 




(91) 
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Figure 11: Rapidity distributions obtained by changing 1/2 < kp < 2 where kp = yU/j/yUF- 
Here we choose Q = Mz = fip- The evolution is based on Candia using MRST grid 
input with fiQ = 1.25 GeV^. As before Vrap has been used for the calculation of the 
hard scattering 
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Figure 12: Plot of the rapidity distributions at LO, NLO and NNLO for Alekhin's model 
and MRST. Shown are also the bands due to the variation of the iip scale, and the errors 
on the cross sections at the corresponding orders. 
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the Alekhin set these equal - for Q = 50 GeV - 2.8% of the best fit value, raising to almost 
3% at Mz and decreasing to 1.7% at 200 GeV. For the MRST set we have determined the 
error on the NLO cross sections in Tab. [H They are about 2.8% of the best fit value at 
Q = 50 GeV, decrease to 1.6% on the peak and decrease moving toward the tail, equating 
1.2% at 200 GeV. 

A more complete view of the role played both by the errors at each perturbative order 
and the corresponding best-fit values can be obtained from Figs. [T3lfT6l In Fig. [13] we 
have zoomed on the region of invariant mass of the lepton pair around 100-102 GeV 
and presented plots of the Alekhin model with the relative errors. The NLO and NNLO 
predictions show overlapping error bands, while the NLO error band for the MRST set 
(Fig. [HI) lays slightly below the Alekhin's result at the corresponding order. We have also 
tried to provide an overall view of the tail of the distribution in Fig. [151 where we show 
the best-fit result at NNLO and the NLO error band. The best-fit value lays outside this 
band. A similar result holds for the MRST set and is shown in Fig. [TJl Finally, in Fig. [15] 
we have zoomed over the region of the resonance, where the best-fit value is shown to lay 
inside the error band. 



8 Estimating the size of the QCD corrections due to 
the evolution 

To estimate the role played by the evolution in determining the full NNLO prediction, 
we show in three tables results for some approximations of the NNLO DY cross sections 
obtained by varying either the hard scattering or the order of the evolved pdf's in the 
factorization formula. These approximations may serve as possible ways to estimate the 
contribution coming from the evolution from that of the hard scatterings, and may provide 
some partial information on their role in the final result. Tables [TT] and [18] show that the 
error made by neglecting the NNLO corrections in the hard scatterings - while keeping 
the entire NNLO evolution - is around 2-3 %, and the correct NNLO cross section is both 
underestimated and overestimated, while a slight bigger error is made if we neglect the 
NNLO corrections to the evolution (4%). The overall decrease of the total cross section 
appears only after the inclusion of the NNLO evolution. In a final table (Tab. [T9j) we 
repeat the trick at NLO, by keeping the hard scatterings at NLO and convoluting with 
the NNLO evolution. Also in this case the errors are around 4% and below in the region 
of Q that we have studied. 

There are some features which are quite evident from this analysis. The first is that 
both at NLO (tab. ^ and at NNLO (tab. ^ the role of the NNLO terms is to reduce 
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the contribution to the cross section. A second piece of information can be extracted by 
comparing all the tables, and extracting the differences Aq = <Jnnlo ® ^nnlo — o'nlo ® 
^NNLO over the entire range of variability of Q and comparing them with the canonical 
NLO cross section, ct^lo'^^nlo- One can easily come to the conclusion that by combining 
the NNLO evolution with the NLO hard scatterings this "improved" NLO cross section 
is much closer to the true NNLO result than the canonical NLO approximation obtained 
using the NLO pdf's. For instance for Q=50 GeV the improved NLO result differs by 
0.3% from the correct determination, while the ordinary NLO prediction differs from it 
by 4%. On the Z peak the improved NLO result differs by 1% respect to the correct 
NNLO prediction, while the standard NLO cross section is 4% away. The pattern is quite 
general. It would be interesting to test the same approach on other NNLO computations 
and check whether on a more general basis the NLO "improved" cross section can be 
used also for other processes as a better estimate of the NNLO result when this is not 
available. We have seen that using these types of approaches, one can estimate the role 
played by the evolution, which in DY dominates over the NNLO corrections to the hard 
scatterings. 

It is clear from the results of these studied that the role played by the NNLO QCD 
corrections in the i^'-factors at NNLO on the Z peak is relevant, corresponding to varia- 
tions that can be reasonably assumed around the few percent level. We recall that with 
10~^ fb of integrated luminosity the statistical error expected on the Z peak is around 
0.05% at the LHC As we have mentioned, suitable choices of the electroweak parameters 
allow to take into account the bulk of the electroweak effects, while the non-factorizable 
contributions are not included in this approach [23] • It is then clear that, given the size 
of the QCD NNLO corrections, we need to worry about these additional effects, which 
are clearly dominant especially if we are interested in having a robust determination of 
all the contributions to this process. Searching for heavy extra Z' is going to be critically 
linked to the correct quantification of these additional corrections |23j . 

9 Conclusions 

We have presented a comparative study of the NNLO predictions for lepton pair pro- 
duction and discussed their robustness. We have presented results concerning i^-factors, 
renormalization/factorization scale dependence and errors on the cross sections - induced 
by errors on the pdf's - following different approaches. For this reason we have put under 
close scrutiny the theory of the logarithmic expansions, which we have shown to give 
results which are compatible with other approaches, and allows to address the issue of 
accuracy in the context of the QCD evolution. Our estimate of the difference between the 
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da/dQ in [pb/GeV] for Alekhin with Q'^^^x^^ VS ^ U TeV 


Q [GeV] 




O'NLO 


(^NNLO 


50 


6.22 ± 0.27 


7.48 ± 0.24 


7 A3 ± 0.21 


60.04 


3.72 ± 0.15 


4.50 ± 0.13 


4.49 ± 0.12 


70.1 


3.30 ± 0.12 


4.03 ± 0.11 


4.05 ± 0.10 


80.1 


6.65 ± 0.24 


8.20 ± 0.24 


8.19 ± 0.23 


90.19 


253 ± 8 


313 ± 9 


309 ± 8 


91.19 


415 ± 14 


514 ± 15 


506 ± 15 


120.07 


0.80 ± 0.02 


0.99 ± 0.02 


0.96 ± 0.03 


146.1 


0.225 ± 0.006 


0.277 ± 0.007 


0.269 ± 0.007 


172.1 


0.097 ± 0.002 


0.119 ± 0.003 


0.117 ± 0.003 


200 


0.047 ± 0.001 


0.058 ± 0.001 


0.058 ± 0.001 



Table 15: Cross sections derived from the best fits for the 3 orders with their errors for 
the set by Alekhin. 



da/dQ in [pb/GeV] for MRST with = ^1 = y/S = U TeV 


Q [GeV] 


cnlo 


50 


6.77 ± 0.19 


60.04 


4.13 ± 0.10 


70.1 


3.79 ± 0.08 


80.1 


7.90 ± 0.14 


90.19 


305 ± 5 


91.19 


499 ± 8 


120.1 


0.952 ± 0.014 


146.1 


0.264 ± 0.003 


172.1 


0.113 ± 0.001 


200 


0.0556 ± 0.0007 



Table 16: Cross sections derived from the best fits at NLO with the errors for the MRST 
set. 
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Figure 13: Errors of the pdf's on the cross sections at LHC. Zoom in the region of 100 
GeV. 
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Figure 14: NNLO cross section for Alekhin with the respective pdf 's errors at NLO in the 
100 GeV region, with Q = fip = fJ-R and \/S = 14 TeV 
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different approaches is slightly above the level of 1%. The factors found using these 
new methods appears to be slightly larger than those coming from the MRST and Alekhin 
evolved parton distributions, but compatible with them, given the actual errors on the 
pdf 's. Clearly, with the advent of the LHC, these analysis should be rendered even more 
accurate, especially at large values of the mass distributions, where a detailed analysis of 
the electroweak effects should be included. This is particularly important in the search of 
extra gauge interactions using this channel. On the Z resonance these effects are smaller, 
at the percent level, but are important for calibration and partonometry. These and other 
related issues will be left for future work. 



Note Added 

The extended analysis presented in this work has been performed within the 2006 
Monte Carlo workshop held in Frascati under the sponsorship of INFN of Italy. Detailed 
tables/plots are provided only for reference in this version for the arXiv since they may 
be of practical use. Candia, CANDlA^y and their interface with Vrap will be released 
and described in forthcoming work. 
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Figure 16: NNLO cross section for Alekhin and MRST with the respective pdf 's errors 
at NLO in the Z peak region. Zoom in the region of 90 GeV. 
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Candia evolution with MRST input, /x^ = 1.25 GeV^ 


Q [GcVl 

^ L J 


^NNLO ® ^NNLO 


^NLO ^ ^NNLO 


Sa 








3^195-1 0-3 


fin n4fiQ 

\J\J •\Jrt\Jc/ 


QQQ7 . 1 n+0 


3 Q8fi4 • 1 0+'' 


3 3040 • 1 fl-^ 




Q fiQfiO . 1 0+*^ 


3 fifi83 • 1 0+*^ 


7 5645 • 10-^ 


80.1407 


7 fi755 • 1 0+" 


7 fifi3Q • 10+" 


1 5087- 10"^ 


on iS7fi 


9 9395 • 1 0+2 


Qfi7fi . 1 n+2 


1 1 988 - 1 0-2 


91.1876 


4 8006 • 10+2 


4 8644 - 10+2 


1 3293 - 10-2 


92 1876 


2 9179 • 10+2 


2 9604 • 10+2 


1 4556 • 10-2 


120.0701 


9.0411 • 10-1 


9.3152 • 10-1 


3.0318-10-2 


146.0938 


2.5267- 10-1 


2.5981 • 10-1 


2.8222 • 10-2 


172.1175 


1.0938-10-1 


1.1179-10-1 


2.2061 - 10-2 


200.0000 


5.4431 • 10-2 


5.5145-10-2 


1.3118- 10-2 



Tabic 17: (Jnnlo ® ^nnlo vs ct^lo ® ^nnlo in [pb/GeV]. Comparison between NNLO 
and NLO cross sections obtained by the convolution of NNLO pdf's. Sa is defined as 

WnNLO <8) ^NNLO — (7 NLO ® ^NNLo\/'^NNLO ^NNLO- 



Candia evolution with MRST input, = 1.25 GeV^ 


Q [GcA'] 






Sa 


50.0000 


6.4935 • 10+° 


6.7853 • 10+" 


4.4938 • 10-2 


60.0469 


3.9997 - 10+" 


4.1805 - 10+" 


4.5206- 10-2 


70.0938 


3.6962 - 10+" 


3.8571 - 10+" 


4.3521 - 10-2 


80.1407 


7.6755 - 10+" 


7.9669 - 10+" 


3.7967-10-2 


90.1876 


2.9325 • 10+2 


3.0219 • 10+2 


3.0498 • 10-2 


91.1876 


4.8006 • 10+2 


4.9437 • 10+2 


2.9809 • 10-2 


92.1876 


2.9179 • 10+2 


3.0029 • 10+2 


2.9141 • 10-2 


120.0701 


9.0411 - 10-1 


9.2025 - 10-1 


1.7858-10-2 


146.0938 


2.5267- 10-1 


2.5593-10-1 


1.2890- 10-2 


172.1175 


1.0938-10-1 


1.1039-10-1 


9.2157-10-3 


200.0000 


5.4431 • 10-2 


5.4781 • 10-2 


6.4302 • 10-3 



Table 18: (Tnnlo ® ^nnlo vs aNNLO ® ^nlo in [pb/GeV]. Upper bound on the NNLO 
cross sections obtained by the convolution of NNLO and NLO pdf's. 5a is defined as 

WnNLO ® ^NNLO — NLO ® ^NLol/'^NNLO ^NNLO- 
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Candia evolution with MRST input, /il = 1.25 GeV^ 


Q [GeV] 


CTNLO ^NLO 


CTNLO ^NNLO 


5a 


50.0000 


6.8119- 10+" 


6.5164 • 10+" 


4.3376 • 10-2 


60.0469 


4.1552 - 10+" 


3.9864- 10+" 


4.0623 - 10-2 


70.0938 


3.8110- 10+" 


3.6683 - 10+" 


3.7465 - 10-2 


80.1407 


7.9371 • 10+" 


7.6639 - 10+" 


3.4420 - 10-2 


90.1876 


3.0657 • 10+2 


2.9676 • 10+2 


3.2000 • 10-2 


91.1876 


5.0242 • 10+2 


4.8644 • 10+2 


3.1790 • 10-2 


92.1876 


3.0569-10+2 


2.9604- 10+2 


3.1584-10-2 


120.0701 


9.5677- 10-1 


9.3152- 10-1 


2.6396-10-2 


146.0938 


2.6562- 10-1 


2.5981 - 10-1 


2.1896-10-2 


172.1175 


1.1382- 10-1 


1.1179-10-1 


1.7801 - 10-2 


200.0000 


5.5940 • 10-2 


5.5145 • 10-2 


1.4212 • 10-2 



Table 19: gnlo ® ^nlo vs aNLO <8) ^nnlo in [pb/GeV]. Lower bound on the NLO 
cross sections obtained by the convolution of NLO and NNLO pdf's. 5a is defined as 

Ic^ATLO ® ^NLO — CFnLO ® ^NNLo\/(^NLO ® ^NLO- 



ctsiQ'^) for MRST at NNLO. Brute force vs Aqcd parameterization. 


Q [GeV] 








50.0000 


0.1268 


0.1251 


1.3702 


83.4897 


0.1170 


0.1156 


1.2570 


90.7209 


0.1156 


0.1142 


1.2405 


9L1876 


0.1155 


0.1141 


1.2395 


92.0543 


0.1153 


0.1139 


1.2377 


107.0583 


0.1128 


0.1115 


1.2090 


125.6466 


0.1103 


0.1090 


1.1801 


144.2350 


0.1082 


0.1069 


1.1564 


162.8233 


0.1064 


0.1052 


1.1363 


181.4117 


0.1049 


0.1037 


1.1191 


200.0000 


0.1035 


0.1024 


1.1040 



Table 20: NNLO running of the couphng determined using the brute force solution of the 
renormalization group equations versus the asymptotic expansions in terms Aqcd- The 
percentage differences are normalized respect to ahmte- 
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xuv{x) Candia evolution at NLO, Les Houches input, Nf = Q = /j,f = /J'R = 100 GeV 


X 


/ N Candia 

Til ^. ( T 


K = 1 


K = 2 


K = 3 


K = 4 


^ = 5 


K ~ 6 


le — Do 


z.Do/o • lU 


CI 01 1 n— 3 
z.oizi • lU 


/I con 1 n^3 


A coc 1 n^3 
Z.40Z0 - iU 


O /I oqq 1 n^3 
Z.4000 - iU 


0/1QQ1 in— 3 
z.4ool • iU 




n nnm 
U.UUUi 


1 oc/1 /I 1 n— 2 


1 O/l 7/1 1 n — 2 
i.Z4 1 4 • iU 


1 ooon 1 n— 2 
i.zzyu • iU 


1 oqn7 1 n— 2 


1 oqrif^ 1 n— 2 
i.zoUD - iU 


1 OQH/^ 1 n— 2 
i.zoUO • iU 


i.ZoUO • iU 


U.UUi 


r '70Q'7 1 n— 2 

0. (ao ( • iU 


p; 7QOQ 1 n— 2 
0. lo\)6 ■ iU 


71 ^1 1 n— 2 
0. / ioi • iU 


r fTo^n 1 n— 2 
0. / zoU • iU 


c ^f^Ac 1 n— 2 
0. / z4D • iU 


0. / z4o • iU 


0. / z4o • iU 


U.Ui 






o Qi S2Q 1 n— 1 


o Q91 Q 1 n— 1 


9 Q9no 1 n— 1 
z.ozuy • iU 




z.ozuy • iU 


U.i 


0.0400 • iU 


r crccc 1 A— 1 
O.OOOD • iU 


0.0 ( oo • iU 


r CfT-l o 1 n^i 
o.O ( iz • iU 


K c:7i R 1 n— 1 
0.0/ ID • lU 


0.0/ 10 • iU 


0.0/ iO • iU 


n 


/I 01 nf^ 1 n— 1 


4.o4y4 • lU 


/I QQQA in — 1 


/I Q7Q/1 in^i 
4.0 I o4 - iU 


A Q700 1 n— 1 
4.0 ( yz - iU 




A in— 1 
4.0 / yi • iu 


yj.o 


o.ooyo • lU 


q A cnq i 1 


q /)Qq7 1 fi— 1 


q /i7QP^ 1 n^i 
0.4/80 - iU 


q A 70 q 1 n^i 

0.4 ( yo - iu 


Q A 7no 1 n— 1 
0.4 r yz • iU 


o.4/yz • iu 


U.4 


oQft/i 1 n— 1 

Z.ZoU4 • iU 


1/1 7n 1 n— 1 


1 70O 1 n— 1 
z. i < zy • lU 


1 AQO 1 n— 1 
z. iooy - lu 


1 RC\K 1 n— 1 
z.ioyo - iu 


z.ioy4 • iu 


Z.iDy4 • iU 


u.o 


1 9971 in~i 


i.lDDl • iU 


1 18*^(1 in— 1 


1 1 sn'^ in— 1 


1 1 sns in— 1 


1 1007 in~i 

i.ioU / • iU 


l.ioU / • iU 


0.6 


5.6866 • 10-2 


5.3292 • 10-2 


5.4212 • 10-2 


5.4067 • 10-2 


5.4090 • 10-2 


5.4086 • 10-2 


5.4087 • 10-2 


0.7 


2.0429-10-2 


1.8840-10-2 


1.9232 • 10-2 


1.9169-10-2 


1.9179-10-2 


1.9178- 10-2 


1.9178 - 10-2 


0.8 


4.6754 • 10-3 


4.2230 • 10-3 


4.3329-10-3 


4.3152 - 10-3 


4.3180- 10-3 


4.3175 - 10-3 


4.3176 - 10-3 


0.9 


3.6098 • lO-'i 


3.1538- lO-'i 


3.2674- 10-1 


3.2486 - 10-1 


3.2516- 10-1 


3.2511- lO-'i 


3.2512-10-'! 



Table 21: NLO valence distribution of the up quark with the Les Houches benchmark 
model 



xuy{x) Candia evolution at NNLO, Les Houches input, N f — A, Q — jxp — — 100 GeV 


X 


/ N Candia 


K = 2 


= 3 


K = 4 


K = 5 


K = & 


K = l 


le-05 


3.0260 • 10-3 


2.9464 • 10-3 


2.8972 • 10-3 


2.8928 • 10-3 


2.8945 • 10-3 


2.8944 • 10-3 


2.8944 • 10-3 


0.0001 


1.3656 • 10-2 


1.3320 • 10-2 


1.3179-10-2 


1.3177-10-2 


1.3181-10-2 


1.3181-10-2 


1.3181 • 10-2 


0.001 


5.9360 - 10-2 


5.8657-10-2 


5.8425 - 10-2 


5.8448 - 10-2 


5.8451 - 10-2 


5.8450-10-2 


5.8450-10-2 


0.01 


2.3139-10-1 


2.3254-10-1 


2.3248-10-1 


2.3253-10-1 


2.3252- 10-1 


2.3252- 10-1 


2.3252 - 10-1 


0.1 


5.5125-10-1 


5.5406-10-1 


5.5451 - 10-1 


5.5446-10-1 


5.5446 - 10-1 


5.5446- 10-1 


5.5446-10-1 


0.2 


4.8672 • 10-1 


4.8416 • 10-1 


4.8463 • 10-1 


4.8454 • 10-1 


4.8454 • 10-1 


4.8454 • 10-1 


4.8454 • 10-1 


0.3 


3.5017-10-1 


3.4470 • 10-1 


3.4507 • 10-1 


3.4497 • 10-1 


3.4498 • 10-1 


3.4498 • 10-1 


3.4498 • 10-1 


0.4 


2.2030 • 10-1 


2.1460-10-1 


2.1486-10-1 


2.1478 • 10-1 


2.1479-10-1 


2.1479- 10-1 


2.1479-10-1 


0.5 


1.2099-10-1 


1.1660-10-1 


1.1676-10-1 


1.1671 - 10-1 


1.1671 - 10-1 


1.1671- 10-1 


1.1671 - 10-1 


0.6 


5.5957-10-2 


5.3309 - 10-2 


5.3392 - 10-2 


5.3364 - 10-2 


5.3366 - 10-2 


5.3367-10-2 


5.3367-10-2 


0.7 


2.0052-10-2 


1.8854-10-2 


1.8888-10-2 


1.8877-10-2 


1.8877-10-2 


1.8878-10-2 


1.8878-10-2 


0.8 


4.5726 • 10-3 


4.2288 • 10-3 


4.2387 • 10-3 


4.2352 • 10-3 


4.2355 • 10-3 


4.2356 • 10-3 


4.2356 • 10-3 


0.9 


3.5079 • 10-^ 


3.1629 • 10-^ 


3.1736 • 10-^ 


3.1698 • 10-^ 


3.1701 • 10-^ 


3.1702-10-^ 


3.1702-10-^ 



Table 22: NNLO Valence distribution for the up quark with the Les Houches benchmark 
model 
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xg{x) Candia evolution at NLO, Les Houches input, Nf = 4, Q = fip = fin = 100 GeV 


.1; 


,,.)('AXDIA 


= 1 


= 2 


= 3 


/, = 4 




H = (i 


le-05 


2.3578 • 10+2 


2.2829 • 10+2 


2.3693 • 10+2 


2.3560 • 10+2 


2.3580 • 10+2 


2.3577-10+2 


2.3578 • 10+2 


0.0001 


9.3027- 10+1 


9.1414-10+1 


9.3237-10+1 


9.2998-10+1 


9.3031 - 10+1 


9.3026 - 10+1 


9.3027-10+1 


0.001 


3.1540- 10+1 


3.1320-10+1 


3.1562 - 10+1 


3.1537-10+1 


3.1540- 10+1 


3.1540- 10+1 


3.1540-10+1 


0.01 


8.1120- 10+° 


8.1098-10+° 


8.1116-10+° 


8.1122-10+° 


8.1120- 10+° 


8.1120- 10+° 


8.1120-10+° 


0.1 


8.9872 • IQ-i 


9.0284 • IQ-i 


8.9826 • 10-1 


8.9877-10-1 


8.9871 • 10-1 


8.9872 • 10-1 


8.9872 • 10-1 


0.2 


2.5540 • IQ-i 


2.5695 • 10-1 


2.5523 • 10-1 


2.5542 • 10-1 


2.5539-10-1 


2.5540 • 10-1 


2.5540 • 10-1 


0.3 


8.3414- 10-2 


8.4005 • 10-2 


8.3353 - 10-2 


8.3422 • 10-2 


8.3414-10-2 


8.3415-10-2 


8.3414- 10-2 


0.4 


2.7017- 10-2 


2.7232 - 10-2 


2.6995 - 10-2 


2.7020 - 10-2 


2.7017-10-2 


2.7017- 10-2 


2.7017- 10-2 


0.5 


8.0411 - 10--'' 


8.1131 - 10--'' 


8.0338 - 10-3 


8.0420 - 10-3 


8.0410-10-3 


8.0411 - 10-3 


8.0411 - 10-3 


0.6 


2.0343- 10-3 


2.0552- 10-3 


2.0321 - 10-3 


2.0345 - 10-3 


2.0342-10-3 


2.0343 - 10-3 


2.0343- 10-3 


0.7 


3.8654 • 10-4 


3.9132- 10-4 


3.8603 • 10-4 


3.8660 • 10-1 


3.8653 • 10-4 


3.8654 • 10-4 


3.8654 • 10-4 


0.8 


4.1712-10-5 


4.2399 • 10-5 


4.1634 • 10-5 


4.1722 • 10-5 


4.1710-10-5 


4.1712-10-5 


4.1712-10-5 


0.9 


1.8310- 10-<^ 


1.8598 • 10-*^ 


1.8272 • 10-<^ 


1.8315- 10-*^ 


1.8309 • 10-6 


1.8310-10-6 


1.8310- 10-6 



Table 23 



NLO gluon density in the Les Houches model 



xg{x) Candia evolution at NNLO, Les Houches input, Nf = 4, Q = fip ^ fi^ = 100 GeV 



X 


/ ^(JANDIA 

.J\ /asymp 


K = 2 


K = 3 


K = 4 


K = 5 


K = 6 


K = 7 


le-05 


2.2328- 10+2 


2.2090 - 10+2 


2.2351 - 10+2 


2.2330 - 10+2 


2.2328 - 10+2 


2.2328- 10+2 


2.2328 - 10+2 


0.0001 


9.0763-10+1 


9.0383 - 10+1 


9.0795-10+1 


9.0765-10+1 


9.0762- 10+1 


9.0763-10+1 


9.0763-10+1 


0.001 


3.1371 • 10+1 


3.1336 • 10+1 


3.1372 • 10+1 


3.1371 • 10+1 


3.1371 • 10+1 


3.1371 • 10+1 


3.1371 • 10+1 


0.01 


8.1407- 10+° 


8.1433 • 10+° 


8.1403 • 10+° 


8.1407-10+° 


8.1407-10+° 


8.1407-10+° 


8.1407 • 10+° 


0.1 


9.0545-10-1 


9.0662- 10-1 


9.0538-10-1 


9.0545 - 10-1 


9.0545 - 10-1 


9.0545-10-1 


9.0545-10-1 


0.2 


2.5753-10-1 


2.5793-10-1 


2.5750-10-1 


2.5752 - 10-1 


2.5753-10-1 


2.5753-10-1 


2.5753-10-1 


0.3 


8.4120- 10-2 


8.4266- 10-2 


8.4112-10-2 


8.4119-10-2 


8.4120-10-2 


8.4120- 10-2 


8.4120- 10-2 


0.4 


2.7238 • 10-2 


2.7288 • 10-2 


2.7235 • 10-2 


2.7238 • 10-2 


2.7238 • 10-2 


2.7238 • 10-2 


2.7238 • 10-2 


0.5 


8.1019- 10-3 


8.1176-10-3 


8.1009 • 10-3 


8.1018-10-3 


8.1019-10-3 


8.1019-10-3 


8.1019- 10-3 


0.6 


2.0476 • 10-3 


2.0518-10-3 


2.0473 • 10-3 


2.0476 • 10-3 


2.0477-10-3 


2.0476 • 10-3 


2.0476 • 10-3 


0.7 


3.8845-10-4 


3.8929 - 10-4 


3.8837-10-4 


3.8845 - 10-4 


3.8845 - 10-4 


3.8845 - 10-4 


3.8845 - 10-4 


0.8 


4.1738- 10-5 


4.1842 - 10-5 


4.1724-10-5 


4.1738-10-5 


4.1738-10-5 


4.1738- 10-5 


4.1738-10-5 


0.9 


1.8861- 10-° 


1.8899-10-° 


1.8853-10-° 


1.8862 - 10-° 


1.8861 - 10-° 


1.8861 - 10-° 


1.8861- 10-° 



Table 24: NNLO gluon density in the Les Houches model 
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xuy{x) Candia evolution at NLO, MRST input, /Uo = 1 GeV, Q = hf = Hr = 100 GeV 


X 


1 -.(Jandia 


K=l 


K = 2 


K = 3 


K = 4 


K = 5 


K = 6 


le - 0") 


l.;-!874 ■ 1()-- 


1.4223 ■ 1()-- 


1.4015 ■ 1U-- 


1.4049 ■ l()-2 


1.404:-! ■ i()-2 


1.4044 ■ l()-2 


1.4044 ■ 10-2 


0.0001 


2.9241 • 10-2 


3.0230-10-2 


2.9654-10-2 


2.9753-10-2 


2.9734- 10-2 


2.9738-10-2 


2.9737-10-2 


0.001 


7.2512 • 10-2 


7.5810-10-2 


7.4246 - 10-2 


7.4530 - 10-2 


7.4472 - 10-2 


7.4485 - 10-2 


7.4482 - 10-2 


0.01 


2.1394-10-1 


2.2294-10-1 


2.2014-10-1 


2.2068-10-1 


2.2057- 10-1 


2.2059- 10-1 


2.2059 - 10-1 


0.1 


5.1946-10-1 


5.1264 - 10-1 


5.1673- 10-1 


5.1598- 10-1 


5.1613- 10-1 


5.1610- 10-1 


5.1611 - 10-1 


0.2 


4.7098 • 10-1 


4.4773 • 10-1 


4.5514-10-1 


4.5372 • 10-1 


4.5402 • 10-1 


4.5395 • 10-1 


4.5397-10-1 


0.3 


3.3696 • 10-1 


3.1056 • 10-1 


3.1763-10-1 


3.1627-10-1 


3.1655-10-1 


3.1649-10-1 


3.1650-10-1 


0.4 


2.0581 - 10-1 


1.8429-10-1 


1.8952 - 10-1 


1.8851 - 10-1 


1.8872-10-1 


1.8868-10-1 


1.8869-10-1 


0.5 


1.0726-10-1 


9.3306-10-2 


9.6490 - 10-2 


9.5870 - 10-2 


9.6001 - 10-2 


9.5972-10-2 


9.5979 - 10-2 


0.6 


4.5848 - 10-2 


3.8656-10-2 


4.0225-10-2 


3.9918-10-2 


3.9983 - 10-2 


3.9968-10-2 


3.9972 - 10-2 


0.7 


1.4636-10-2 


1.1901 - 10-2 


1.2479-10-2 


1.2365-10-2 


1.2390-10-2 


1.2384 - 10-2 


1.2385 - 10-2 


0.8 


2.7935 ■ 10--' 


2.1GG9 ■ 10--' 


2.29G8 ■ 10-"' 


2.271U ■ 10-"' 


2.27()3 ■ K)-' 


2.2753 ■ 10--' 


2.2756 ■ 10-3 


0.9 


1.5488 • 10-^ 


1.1138-10-^ 


1.2032 • 10-4 


1.1852-10-4 


1.1891-10-4 


1.1882- 10-4 


1.1884-10-4 



Table 25: NLO distribution of the quark with MRST input evolved with Candia 



xu^{x) Candia evolution at NNLO, MRST input, /xq = 1 GeV, Q = fip = jjm = im GeV 


X 


/ N Candia 

xu^yx}^^,^,^^ 


K = 2 


K = 3 


K = 4 


K = 5 


K = 6 


K = 7 


le-05 


1.2678-10-2 


1.3052 - 10-2 


1.2939-10-2 


1.2940-10-2 


1.2945- 10-2 


1.2944- 10-2 


1.2944-10-2 


0.0001 


2.9047-10-2 


2.9858-10-2 


2.9616-10-2 


2.9627-10-2 


2.9636 - 10-2 


2.9634- 10-2 


2.9634 - 10-2 


0.001 


7.3813-10-2 


7.5701 - 10-2 


7.5337-10-2 


7.5391 - 10-2 


7.5399 - 10-2 


7.5395 - 10-2 


7.5395 - 10-2 


0.01 


2.1492 ■ 10-1 


2.2015 ■ 10-1 


2.1997 ■ lU-i 


2.2009 ■ 1()-'- 


2.2008 ■ l()-i 


2.2007- l()-i 


2.2007 ■ 10-1 


0.1 


5.2810-10-1 


5.2525-10-1 


5.2612-10-1 


5.2599-10-1 


5.2597-10-1 


5.2598-10-1 


5.2598-10-1 


0.2 


4.8293-10-1 


4.6760 - 10-1 


4.6869-10-1 


4.6841 - 10-1 


4.6841 - 10-1 


4.6842- 10-1 


4.6842 - 10-1 


0.3 


3.4556-10-1 


3.2689-10-1 


3.2781 - 10-1 


3.2753-10-1 


3.2753-10-1 


3.2755-10-1 


3.2755 - 10-1 


0.4 


2.0973-10-1 


1.9409-10-1 


1.9473-10-1 


1.9451 - 10-1 


1.9452 - 10-1 


1.9453- 10-1 


1.9453-10-1 


0.5 


1.0794 • 10-1 


9.7714 • 10-2 


9.8084 • 10-2 


9.7952 • 10-2 


9.7958 • 10-2 


9.7964 • 10-2 


9.7963 • 10-2 


0.6 


4.5227 • 10-2 


3.9989 • 10-2 


4.0165 • 10-2 


4.0101 • 10-2 


4.0104-10-2 


4.0107- 10-2 


4.0106-10-2 


0.7 


1.4011 - 10-2 


1.2057-10-2 


1.2121 - 10-2 


1.2097-10-2 


1.2098-10-2 


1.2099- 10-2 


1.2099-10-2 


0.8 


2.5508 - 10-3 


2.1205 - 10-3 


2.1343-10-3 


2.1291 - 10-3 


2.1294- 10-3 


2.1297- 10-3 


2.1296 - 10-3 


0.9 


1.2957-10-4 


1.0200-10-4 


1.0290-10-4 


1.0257- 10-4 


1.0258 - 10-4 


1.0260- 10-4 


1.0260-10-4 



Table 26: NNLO Uy quark distribution evolved with Candia using the MRST input. 
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xg{x) Candia evolution at NLO, MRST input, //o = 1 GcV, Q = /i_F = I-ir = 100 GeV 


.1; 


(,.)C'AXDIA 
■'-.V V / as y nip 


= 1 


= 2 


= 3 


/, = 4 




H = (i 


le-05 


1.9288-10+2 


1.8441 • 10+2 


1.9429 • 10+2 


1.9262 • 10+2 


1.9293-10+2 


1.9287-10+2 


1.9288 • 10+2 


0.0001 


8.1300-10+1 


7.9222 - 10+1 


8.1613-10+1 


8.1245 - 10+1 


8.1310- 10+1 


8.1297- 10+1 


8.1300-10+1 


0.001 


2.9001 - 10+1 


2.8675 - 10+1 


2.9043 - 10+1 


2.8995 - 10+1 


2.9002- 10+1 


2.9001 - 10+1 


2.9001 - 10+1 


0.01 


7.8335-10+° 


7.8282 - 10+° 


7.8328 - 10+° 


7.8339 • 10+° 


7.8334 - 10+° 


7.8336 - 10+° 


7.8335-10+° 


0.1 


9.3962- IQ-i 


9.4520 • IQ-i 


9.3881 • 10-1 


9.3976 • 10-1 


9.3959-10-1 


9.3963-10-1 


9.3962 • 10-1 


0.2 


2.7632 • IQ-i 


2.7832 • IQ-i 


2.7605 • 10-1 


2.7636 • 10-1 


2.7631 • 10-1 


2.7632-10-1 


2.7632 • 10-1 


0.3 


9.1622- 10-2 


9.2347-10-2 


9.1533-10-2 


9.1635-10-2 


9.1619-10-2 


9.1622-10-2 


9.1622- 10-2 


0.4 


2.9629- 10-2 


2.9876-10-2 


2.9600 - 10-2 


2.9633 - 10-2 


2.9628-10-2 


2.9629 - 10-2 


2.9629- 10-2 


0.5 


8.6249- 10--'' 


8.6999-10-3 


8.6170-10-3 


8.6260 - 10-3 


8.6248 - 10-3 


8.6249 - 10-3 


8.6249 - 10-3 


0.6 


2.0714- 10-3 


2.0901- 10-3 


2.0697-10-3 


2.0716-10-3 


2.0714-10-3 


2.0714- 10-3 


2.0714- 10-3 


0.7 


3.5704 • 10-4 


3.6038 • 10-4 


3.5677 • 10-4 


3.5707-10-1 


3.5704 • 10-4 


3.5704 • 10-4 


3.5704- 10-4 


0.8 


3.2860 • IQ-^ 


3.3189-10-5 


3.2839 • 10-5 


3.2861 • lO-'"^ 


3.2861 • 10-5 


3.2860-10-5 


3.2860 • 10-5 


0.9 


6.4223 • IQ-^ 


6.4990 • 10-^ 


6.4192 • 10-^ 


6.4218 • 10-^ 


6.4225 • 10-^ 


6.4222 • 10-^ 


6.4223 • 10-^ 



Table 27: NLO gluon density with MRST input, evolved with Candia 



xg{x) Candia evolution at NNLO, MRST input, = 1 GcV, Q = Mfl = 100 GeV 


X 


„„^„^CANDIA 
•^yy-^ J asymp 


K = 2 


K = 3 


K = 4 


K = 5 


K ~ 6 


K = 7 


le-05 


1.6068- 10+2 


1.5727-10+2 


1.6089 - 10+2 


1.6075 - 10+2 


1.6063- 10+2 


1.6066 - 10+2 


1.6066-10+2 


0.0001 


7.1188- 10+1 


7.0542 - 10+1 


7.1229-10+1 


7.1197-10+1 


7.1176- 10+1 


7.1181- 10+1 


7.1181-10+1 


0.001 


2.6591 • 10+^ 


2.6526 • 10+1 


2.6596 • 10+1 


2.6591 • 10+1 


2.6589 • 10+1 


2.6590 • 10+1 


2.6590 • 10+1 


0.01 


7.5377 • 10+° 


7.5426 • 10+° 


7.5373 • 10+° 


7.5374 • 10+° 


7.5377-10+° 


7.5376 • 10+° 


7.5376 • 10+° 


0.1 


9.8748-10-1 


9.8964-10-1 


9.8738-10-1 


9.8748 - 10-1 


9.8755 - 10-1 


9.8753-10-1 


9.8753-10-1 


0.2 


3.0191 - 10-1 


3.0265-10-1 


3.0188-10-1 


3.0192 - 10-1 


3.0194-10-1 


3.0193-10-1 


3.0193-10-1 


0.3 


1.0210- 10-1 


1.0235-10-1 


1.0210-10-1 


1.0211-10-1 


1.0211-10-1 


1.0211- 10-1 


1.0211- 10-1 


0.4 


3.3271 • 10-2 


3.3345 • 10-2 


3.3270 • 10-2 


3.3274 • 10-2 


3.3275 • 10-2 


3.3275 • 10-2 


3.3275 • 10-2 


0.5 


9.6593 • 10-3 


9.6776 • 10-3 


9.6592 • 10-3 


9.6606 • 10-3 


9.6607-10-3 


9.6607-10-3 


9.6607- 10-3 


0.6 


2.2882 • 10-3 


2.2915 • 10-3 


2.2882 • 10-3 


2.2885 • 10-3 


2.2885 • 10-3 


2.2885 • 10-3 


2.2885 • 10-3 


0.7 


3.8337-10-4 


3.8368 - 10-4 


3.8339 - 10-4 


3.8345 - 10-4 


3.8343 - 10-4 


3.8343 - 10-4 


3.8343 - 10-4 


0.8 


3.3492- 10-5 


3.3487-10-5 


3.3495 - 10-5 


3.3501 - 10-5 


3.3497-10-5 


3.3498 - 10-5 


3.3498 - 10-5 


0.9 


5.8939- lO-'^ 


5.8839 - lO-'^ 


5.8951 - lO-'^ 


5.8960-10-^ 


5.8947-10-^ 


5.8950-10-^ 


5.8950- lO-'^ 



Table 28: NNLO gluon density evolved with Candia using the MRST input. 
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da^^'^/dQ [pb/GcV] with MRST input, fij = 1.25 GeV^, Candia evolution, ^ = Q, ^/S = U TeV 


Q 


o'NLO asyni. 


K = 1 


K = 2 


K = 3 


K = 4 


K = 5 


K = 6 




6 8121 • in+° 


7 004"^ • 1 0+" 


6 7706 • 1 n+° 

(J.I i \J\J J-\J 


6 8020 • 1 0+° 


6 7972 • 1 n+° 


fi 7980 • 1 0+° 


fi 7978 • 1 0+° 


6f) 0469 


4 1554 • 10+° 


A 2645 • 10+° 


4.1317 • 10+° 


4 1495 • 10+° 


4 1468 • 10+° 


4.1473 • 10+° 


4.1472 • 10+° 


70 0938 


3 811 2 • 10+" 


3 9073 • 10+" 


3 7916 - 10+" 


3 8071 • 10+" 


3 8047 - 10+" 


3 8051 • 10+" 


3 8051 - 10+" 


80.1407 


7 9374 • 1 0+" 


8 1327 • 10+" 


7 9013 - 10+" 


7 9327 - 10+" 


7 9278 - 10+" 


7 9287 - 10+" 


7 9285 - 10+" 


90.1876 


3.0658 • 10+2 


3.1385 • 10+2 


3.0529 - 10+2 


3.0646 - 10+2 


3.0627 - 10+2 


3.0630 - 10+2 


3.0630 - 10+2 


91.1876 


5.0243 • 10+2 


5.1429-10+2 


5.0033 - 10+2 


5.0223- 10+2 


5.0193-10+2 


5.0198-10+2 


5.0197-10+2 


120.0701 


9.5681 • IQ-i 


9.7592 • 10-1 


9.5264 • 10-1 


9.5582-10-1 


9.5531 • 10-1 


9.5540 • 10-1 


9.5538 • 10-1 


146.0938 


2.6563 • IQ-i 


2.7026 • 10-1 


2.6442 • 10-1 


2.6522 • IQ-i 


2.6509-10-1 


2.6511- 10-1 


2.6511 • 10-1 


172.1175 


1.1382 • 10-1 


1.1558-10-1 


1.1328-10-1 


1.1360-10-1 


1.1355-10-1 


1.1356-10-1 


1.1355-10-1 


200.0000 


5.5942 • 10-2 


5.6694-10-2 


5.5662 • 10-2 


5.5805 • 10-2 


5.5782-10-2 


5.5786-10-2 


5.5785-10-2 



Table 29: Drell-Yan cross section at NLO computed with the MRST parametric input 
(//q = 1-25 GeV^) and the evolution performed using Candia. Shown are the cross 
sections for the truncated solutions and the asymptotic cross section. 



^j^yyu) i^iQ [pi,/G(>V] witli -MRST input. //^ = 1.25, G(-V-' ("axdia c\-olution, = //j. = Q. VS = 14 T(>V 


Q 


(^NNLO asym. 


K = 2 


K = 3 


/t = 4 


K = 5 


K = 6 


K = 7 


50.0000 


6.4940 - 10+" 


6.5052 - 10+" 


6.4758 - 10+" 


6.4807- 10+" 


6.4805 - 10+" 


6.4803 - 10+" 


6.4804 - 10+" 


60.0469 


3.9989 - 10+" 


4.0040 - 10+" 


3.9886 - 10+" 


3.9911 - 10+" 


3.9911 - 10+" 


3.9910- 10+" 


3.9910- 10+" 


70.0938 


3.6948 - 10+° 


3.6995 - 10+" 


3.6868 - 10+" 


3.6888 - 10+" 


3.6888 - 10+" 


3.6887- 10+" 


3.6887- 10+" 


80.1407 


7.6740 • 10+° 


7.6871 • 10+° 


7.6600 • 10+° 


7.6641 • 10+° 


7.6642 • 10+° 


7.6640 • 10+° 


7.6640 • 10+° 


90.1876 


2.9335 • 10+2 


2.9395 • 10+2 


2.9283 • 10+2 


2.9299 • 10+2 


2.9299 • 10+2 


2.9299 • 10+2 


2.9299 • 10+2 


91.1876 


4.8027 • 10+2 


4.8124-10+2 


4.7940 • 10+2 


4.7966 • 10+2 


4.7967 • 10+2 


4.7966 • 10+2 


4.7966 • 10+2 


120.0701 


9.0552-10-1 


9.0678-10-1 


9.0326-10-1 


9.0373-10-1 


9.0376-10-1 


9.0373-10-1 


9.0373-10-1 


146.0938 


2.5318-10-1 


2.5334-10-1 


2.5243- 10-1 


2.5255- 10-1 


2.5256- 10-1 


2.5255 - 10-1 


2.5255 - 10-1 


172.1175 


1.0963-10-1 


1.0963-10-1 


1.0927- 10-1 


1.0931-10-1 


1.0932- 10-1 


1.0931-10-1 


1.0931-10-1 


200.0000 


5.4572 • 10-2 


5.4533 • 10-2 


5.4367-10-2 


5.4388 • 10-2 


5.4390 • 10-2 


5.4389 • 10-2 


5.4389 • 10-2 



Table 30: NNLO Drell-Yan cross section with the MRST initial conditions and the evolu- 
tion performed with Candia. We present the results for the various truncated solutions 
and for the asymptotic one. 
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